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An  Automated  Earth  Temperature  Station 


T.  Kusuda  and  F.  J.  Powell 
National  Bureau  of  Standards 
Washington,  D.  C. 

ABSTRACT 

In  order  to  improve  heat  transfer  information  needed  for  under- 
ground installations,  the  Environmental  Engineering  Section  of  the 
National  Bureau  of  Standards  has  collected  hourly  earth  temperature 
data  for  about  a year  at  its  site  in  Gaithersburg,  Maryland.  Tem- 
peratures were  recorded  at  several  depths  under  five  different  types 
of  earth  surface  conditions;  paved-asphalt  in  color,  paved-painted 
white,  bare  ground,  grass  covered,  and  an  occasionally  mown  area  that 
is  covered  by  naturally  grown  grass.  Hourly  measurements  of  earth 
temperatures  at  selected  depths,  ground  surface  heat  flux,  soil  mois- 
ture content,  ambient  air  temperature  and  relative  humidity,  solar 
radiation,  wind  speed  and  direction,  and  rainfall  were  collected 
using  an  automatic  data  acquisition  system. 

This  report  describes,  in  detail,  the  instrumentation  system  used. 
Included  are  the  characteristics  of  the  transducers,  the  data  logging 
apparatus,  the  data  transmission  scheme  and  the  data  processing  tech- 
nique using  a computer. 

A later  publication  is  planned  to  present  detailed  data  and 
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results  obtained. 


1 . Xnt  roduc  t ion 


The  use  of  an  automatic  data  acquisition  system  is  currently  be- 
coming increasing  advantageous  because  of  the  accessibility  and  efficien- 
cy of  computers  for  processing  a large  amount  of  data.  There  is  a great 
deal  of  discussion  about  automatic  data  acquisition  systems,  but  many 
environmental  scientists  and  engineers  have  not  been  exposed  to  the 
several  technical  problems  associated  with  this  modern  equipment. 
Manufacturer’s  literature  are  usually  presented  in  summary  form  for 
multi-purpose  use  and  sometimes  confuse  those  who  must  plan  specific 
experimental  projects  using  an  automated  data  acquisition  system. 

Modern  electronic  instruments  are  quite  sensitive  to  their  operating 
environment  and  full  time  attention  of  an  experienced  operator  is  usually 
required.  When  a system  is  operating  correctly  it  will  produce  an 
enormous  amount  of  data  but  if  not  properly  planned  for  it  can  over- 
whelm the  experimental  analysts.  Moreover,  data  acquisition  systems 
are  becoming  more  and  more  sophisticated  with  new  models  constantly 
being  introduced.  The  latest  models  can  incorporate  a computer  pro- 
gramed to  control  the  entire  experiment  including  reduction  of  the 
data  collected. 

The  purpose  of  this  report  is  to  describe  in  detail  a simple  and 
relatively  unsophisticated  automatic  data  acquisition  system  which  was 
used  at  NBS  in  an  earth  temperature  experimental  project.  Although  the 
system  is  now  a few  years  old,  its  basic  principles  and  operation  are 
the  same  as  those  of  today's  more  complicated  and  computerized  systems. 


This  report  is  intended  primarily  for  those  considering  the  use 
of  this  type  of  system  and  is  therefore  expository  in  nature  rather 
than  presenting  a concise  summary  or  only  unique  features  of  an  auto- 
matic data  acquisition  system.  Major  emphasis  is  given  to  the  explana- 
tion of  the  several  specific  aspects  of  automatic  data  handling.  A 
brief  statement  of  the  objectives  of  the  NBS  earth  temperature  measure- 
ment project,  to  which  this  particular  data  acquisition  system  was 
addressed,  is  given  in  the  next  section  as  background. 

2 . Background 

Since  July,  1968,  hourly  observations  of  earth  temperature  and 
surface  weather  data  have  been  collected  at  a site  on  the  grounds  of 
NBS  Gaithersburg  campus  (77°  13'  longitude  and  39°  8'  latitude,  ele- 
vation 42  0'  ). 

The  purpose  of  this  data  gathering  effort  was  to  provide  improved 
earth  temperature  information  pertinent  to  the  design  and  evaluation 
of  the  Office  of  Civil  Defense  (OCD),  Department  of  Defense  below 
grade  installations. 
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Previous  studies  conducted  at  the  National  Bureau  of  Standards 


revealed  that  the  heat  absorbing  capacity  of  the  earth  surrounding  an 
occupied  underground  survival  shelter  is  a crucial  factor  for  the  main- 
tenance of  a tolerable  shelter  temperature.  This  was  found  to  be 
especially  true  during  a summer  season  during  a prolonged  and  "sealed- 
up"  condition.  The  heat  absorbing  capacity  of  the  earth  is  in  turn 
governed  by  its  initial  temperature  at  the  time  of  the  shelter  occupancy 
called  here  its  undisturbed  condition. 

An  extensive  literature  search  was  conducted  for  data  on  earth 
temperatures  and  ground  heat  transfer.  This  search  revealed  that: 

1.  The  availability  of  earth  temperature  data  to  depths 
beyond  5 ft  is  extremely  limited.  There  is  a need 
for  reliable  temperature  data  at  greater  depths  for 
heat  transfer  analysis  of  underground  installations. 

2 . The  effect  of  earth  cover  on  temperature  at  several 
depths  was  not  explored  on  a long  time  basis.  It 
was  anticipated  that  there  may  be  an  effective  way 
of  controlling  the  earth  temperature  by  the  proper 
selection  of  earth  surface  covers. 
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3.  Correlations  involving  earth  temperature  changes  with 
surface  weather  characteristics,  such  as  air  temperature, 
humidity,  solar  radiation,  rainfall  and  soil  moisture  are 
rare.  This  is  primarily  due  to  a lack  of  data  recorded 
simultaneously  on  an  hourly  basis  at  any  single  weather 
station  in  the  United  States. 

4.  Earth  temperature  data  currently  published  apply  largely 
to  undisturbed  earth  remotely  located  from  nearby  in- 
stallations. Except  for  a few  cases,  the  effect  of 
nearby  buildings  upon  earth  temperature  has  not  been 
well  explored. 

5.  Data  on  the  effects  of  the  presence  of  ground  water 
upon  the  earth  temperature  was  not  evident  in  the 
literature  studied. 

6.  An  analysis  of  earth  temperature  measurement  error 
for  a practical  probe  (strong  enough  to  withstand 
rough  treatment  during  the  installation  and  trans- 
portation) was  not  revealed. 

7.  Use  of  automated  data  acquisition  and  processing  by 
modern  digital  equipment  for  collection  of  earth 
temperature  data  was  not  reported. 
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8.  As  a result  of  previous  studies  several  unpublished 
comments  indicating  the  presence  of  stray  ground 
electric  currents  affecting  or  overriding  the  actual 
emf  developed  in  thermocouple  circuits  were  revealed. 

The  objective  of  this  investigation  was  to  develop  information 
concerning  the  eight  points  listed  above.  Major  emphasis  was  placed 
on  determining  the  effect  of  earth  cover  on  the  seasonal  variation  of 
temperatures  at  different  depths.  A second,  equally  important,  ob- 
jective was  to  study  of  the  hourly  change  of  near-surface  soil  tem- 
perature as  it  is  affected  by  atmospheric  conditions  above  it. 

A major  reason  for  the  lack  of  this  type  of  information  is  that 
collection  of  a vast  amount  of  data  is  required  over  a period  of  at 
least  one  year  to  cover  all  seasons.  This  entails  a considerable  ex- 
penditure in  terms  of  instrumentation  and  man-hours.  However,  the  in- 
troduction of  computers  and  automatic  data  acquisition  devices  made  it 
feasible  to  undertake  such  experiments  at  less  cost  and  with  fewer 
man-hours  than  would  have  been  required  previously. 
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3.  Description  of  Earth  Temperature  Site 


The  site  designated  for  installing  temperature  probes  was  chosen 
to  be  approximately  500  ft  away  from  the  nearest  building  (Figure  1). 

As  shown  on  a topographic  map  (Figure  2)  the  ground  surface  slopes  down 
toward  west  at  a grade  of  1.25,  providing  good  drainage.  Core  boring 
samples  from  depths  of  1 ft  to  down  to  30  ft  revealed  that  the  soil  is 
of  yellow  silt;  its  average  dry  density  and  moisture  content  being 

3 

about  100  lb/ft  and  15  percent  of  the  dry  weight,  respectively.  The 
moisture  content  of  the  surface  soil  was  higher  than  that  of  the  deep 
soil.  It  was  24  percent  by  weight.  (See  the  soil  test  results  attached.) 

Two  wells  drilled  to  a depth  of  30  ft  in  the  vicinity  of  the  tem- 
perature measurement  site  showed  no  trace  of  ground  water  throughout 
the  test  period.  Five  square  patches  each  50'  x 50'  were  selected  for 
five  different  surface  covering  treatments,  as  indicated  in  Figure  2. 

The  surface  conditions  were:  No.  1,  paved  asphalt  (5"  thick)  with 
natural  asphalt  color;  No.  2,  asphalt  paved  (5"  thick)  and  painted 
white  with  highway  paint  (Federal  Spec.  8010-900-2936);  No.  3,  bare 
ground;  No.  4,  grass  covered  but  periodically  mowed;  and  No.  5,  grass 
naturally  grown.  Each  square  patch  is  separated  from  the  neighboring 
patches  by  a distance  of  50  ft.  A 4"  diameter  x 30’  deep  hole  was  bored 
at  the  geometric  center  of  each  patch.  A temperature  probe  was  installed 
in  each  hole.  The  soil  brought  out  during  the  boring  was  backfilled 
and  compacted  around  the  temperature  probe,  as  close  to  the  condition 
of  the  original  undisturbed  soil  structure  as  was  possible. 
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4.  Earth  Temperature  Probes 


\ 


Each  of  the  five  temperature  probes  contained  twelve  pairs  of  #24 
(AWG)  copper-constantan  thermocouples  assembled  as  shown  in  Figure  3C. 
Each  thermocouple  wire  was  insulated  with  polyvinyl  chloride  plastic 
and  was  manufactured  as  duplex  wire  by  the  Thermoelectric  Company. 

Seven  of  the  twelve  thermocouples  were  centered  in  a water-tight 
polyethylene  tube  of  1"  outside  diameter  and  1/16"  wall  thickness. 

The  thermocouples  were  located  along  each  tube  to  correspond  to  depths 
indicated  in  Figure  3.  Ottawa  sand  was  used  to  fill  the  void  around 
the  thermocouple  in  the  polyethylene  tube.  Five  thermocouples  were 
used  for  depths  less  than  1'  and  branched  out  from  the  main  body  of 
the  probe  to  spread  out  horizontally  around  the  probe  at  depths  as 
illustrated  in  Figure  (3-A). 

In  place,  the  main  vertical  probe  housing  curved  horizontally  at 
the  1'  depth  and  extended  about  51"  to  an  underground  junction  box. 

This  was  done  to  minimize  stem  heat  conduction  error  for  shallow  tem- 
perature measurements.  The  tip  of  each  of  the  five  shallow  thermocouples 
were  protected  from  direct  contact  with  the  soil  by  rubber  sealant  caps. 
This  was  done  to  isolate  the  thermocouple  circuit  from  any  ground  loop 
electric  circuit,  which  would  have  caused  serious  common  mode  problems 
as  described  later.  Protection  of  the  thermocouples  within  the  poly- 
ethylene tubes  and  rubber  sealant  caps  could  cause  error  when  measuring 
earth  temperature  because  such  an  arrangement  prevents  the  thermocouples 


7 


from  direct  contact  with  the  soil.  However,  these  protection  measures 
were  considered  to  be  necessary  to  minimize  not  only  the  effect  of 
ground  loop  current  but  also  to  minimize  mechanical  damage  of  thermo- 
couples during  transportation  and  installation  especially  for  an  ex- 
periment of  this  nature. 

In  order  to  determine  the  measurement  error  due  to  the  protective 
coverings  around  the  thermocouples  a laboratory  investigation  was  made 
on  two  vertical  probes:  one  with  a stainless  steel  tube  and  the  other 
with  a phenolic  tube. 

Figure  4 shows  a sketch  of  the  probes  which  were  separately  in- 
stalled in  a sandbed  as  illustrated  in  Figure  5.  The  sandbed  was  in 
an  18  1/4"  diameter  drum  26"  high  insulated  circumferentially  with  a 
7"  thick  fibrous  glass  blanket  to  simulate  the  adiabatic  boundary  con- 
dition found  in  the  earth.  Figure  5 shows  the  locations  of  the  thermo- 
couples with  those  within  the  protective  probe  matched  exactly  in  depth 
with  those  installed  directly  in  the  surrounding  sand  at  distances  of 
6 1/8"  and  3 1/8"  from  the  center  of  the  probe.  To  simulate  daily 
ground  surface  temperature  fluctuations  two  250  watt  heat  lamps  were 
energized  during  working  hours,  8:30  a.m.  to  5:00  p.m.,  and  were  off 
during  the  remaining  hours  of  the  day.  The  on-off  cycle  was  maintained 
for  a period  of  seven  days  to  attain  a steady  periodic  cycle. 
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Figure  6 shows  sample  results  obtained  from  this  investigation. 
Temperature  data  obtained  from  the  probe  thermocouples  are  compared 
with  data  from  thermocouples  installed  directly  in  the  sand  6 1/8"  away 
from  the  probe.  Figure  6 shows  that  readings  from  the  probe  thermo- 
couples closely  reflect  actual  depth  sand  temperatures  for  depths  more 
than  12".  This  data  is  for  a conductive  stainless  steel  probe  pro- 
truding above  the  sand  surface  and  directly  exposed  to  radiant  heat. 

The  temperature  deviations  between  the  thermocouples  in  the  phenolic 
probe  and  those  in  the  undisturbed  sand  were  found  to  be  even  smaller 
than  those  for  the  stainless  steel  probe  as  given  in  Figure  6A.  It 
is  concluded  that  the  manner  in  which  the  actual  thermocouple  probes 
were  installed,  as  illustrated  in  Figure  3,  yields  very  small  errors 
due  to  the  stem  heat  conduction  and  to  the  thermal  insulating  effect 
of  the  protective  tubing. 

Theoretically,  probe  stem  heat  conduction  error  becomes  smaller 

as  the  thermal  diffusivity  of  the  probe  material  approaches  the  soil 

thermal  diffusivity.  The  typical  soil  thermal  diffusivity  is  0.02 

2 

ft  /hr,  whereas  those  of  stainless  steel,  phenolic,  polyethylene  and 
sand  are  respectively  0.172,  0.004,  0.06  and  0.01. 
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From  the  standpoint  of  thermal  diffusivity  and  for  practical  reasons 
a polyethy lene-sand  filled  probe  was  selected.  A practical  advantage  of 
the  polyethylene  probe  is  that  it  is  not  rigid,  and  thus  can  be  trans- 
ported in  coiled  form  as  compared  with  a long  rigid  rod.  Moreover,  with 
sand  inside  the  tubing  it  becomes  heavy  enough  to  allow  easy  insertion 
into  a vertical  hole  in  the  ground. 

Horizontal  thermocouples  for  shallow  depth  temperature  measurements 
were  installed  at  the  bottom  of  a 1/4"  wide  groove  at  the  depth  indicated 
along  the  rays  shown  in  Figure  3A.  The  narrow  groove  was  then  backfilled 
with  the  surrounding  soil.  In  this  manner  the  earth  was  least  disturbed 
and  the  depths  of  the  thermocouples  were  controlled  very  accurately. 

5.  Junction  Box  and  Underground  Cable 

In  addition  to  thermocouples,  soil  moisture  gauges  and  heat  flow 
transducers  as  described  in  sections  7.6  and  7.3,  respectively,  were 
installed  around  the  thermal  probe  as  shown  in  Figure  3.  The  lead 
wires  from  all  of  the  transducers  for  each  patch  were  brought  to  an 
underground  junction  box.  The  box  was  12"  long,  8"  wide  and  6"  high 
and  was  made  of  galvanized  cast  iron.  This  box  contained  all  waterproof 
connections  to  the  underground  cables  leading  to  the  instrument  bunker. 
Figure  3C  shows  a schematic  diagram  of  junction  box  terminals  for  the 
various  transducers.  Each  junction  box  was  sealed  with  a rubber  gasketed 
lid  and  inlets  to  the  box  were  well  sealed  with  a silicone  rubber  com- 
pound (Dow  Corning  732  RTV  Silastic). 
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Each  underground  cable  contained  21  individually  insulated,  labelled, 
solid  conductors  all  #24  A.W.G.;  16  were  copper  and  5 were  constantan. 

All  wires  were  enclosed  in  a thick  polyvinyl-chloride  cover  to  form  a 
cable  considered  appropriate  for  an  underground  installation.  The  manu- 
facturer of  the  multi-conductor  cable  was  the  Thermoelectric  Company.  A 
cross  section  of  the  underground  cable  is  shown  in  Figure  3C.  The  wires 
were  arranged  for  easy  identification. 

Patch  #3,  bare  ground  surface,  was  provided  with  two  sets  of 
junction  boxes  and  underground  cables,  because  it  contained  7 more 
thermocouples  and  three  more  heat  flow  transducers  than  the  other  four 
patches.  The  extra  thermocouples  were  purposely  installed  unprotected 
and  were  placed  outside  the  polyethylene  probe  to  allow  a check  on  the 
magnitude  of  error  introduced  to  unprotected  thermocouples.  The  standard 
heat  flow  transducers  for  all  patches  were  installed  at  the  4"  depth 
except  at  Patch  #3  as  discussed  in  section  7.3. 

Figure  3D  shows  the  manner  in  which  the  underground  cables  were 
connected  to  the  instrument  bunker  via  24”  deep  ditches. 
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6.  Underground  Instrument  Bunker 


A special  underground  instrument  bunker  was  constructed  to  house 
the  temperature  measuring  station.  Construction  features  and  the  in- 
terior facilities  of  the  bunker  are  schematically  in  Figure  7.  All  but 
the  teletype  unit  of  the  automatic  data  acquisition  system  were  housed 
within  this  bunker.  Some  of  these  instruments  were  temperature  sensitive 
and  the  bunker  was  air  conditioned  during  the  summer  season  to  keep  the 
temperature  below  90  °F  and  relative  humidity  below  507o.  A 1/4  ton 
nominal  capacity  condensing  unit  was  installed  above  the  bunker  with  a 
ceiling  mounted  cooling  coil  inside  the  bunker.  About  one  gallon  per 
week  of  condensate  was  collected  during  a summer  week.  The  bunker 
entry  was  a vertical  hatch  covered  by  a lightweight  manhole  cover. 

A telephone  system  between  the  instrument  bunker  and  a laboratory 
module  at  NBS  was  used  to  transmit  data  from  the  bunker  at  a teletype 
unit  in  the  laboratory.  This  telephone  line  was  also  suitable  for 
voice  communication. 

The  exterior  surfaces  of  the  bunker  walls  were  well  waterproofed 
but  waterproofing  around  the  incoming  bundle  of  six  underground  cables 
was  especially  difficult  to  obtain.  Initially  rain  water  found  its  way 
into  the  bunker  through  this  opening  in  the  wall. 
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Water  leakage  after  the  first  heavy  rain  flooded  the  bunker  floor 
to  a depth  of  12",  causing  considerable  damage  to  several  critical  in- 
struments, such  as  the  data  transmission  sets,  thermocouple  reference 
junctions,  and  wind  data  signal  conditioners.  Several  transformers  in 
other  instruments  were  also  damaged  apparently  due  to  the  effects  of 
this  electrical  storm. 

Thunderstorms  were  also  a cause  of  data  transmission  failure.  A 
lightning  rod  was  later  installed  over  the  instrument  bunker  atop  the 
wind  tower  and  this  helped  decrease  the  number  of  data  transmission 
failures  during  thunderstorms.  Whenever  data  transmission  fails,  the 
system  must  be  manually  restarted. 

7.  Transducer  Systems 

Table  1 shows  data  channel  identification  numbers  to  each  trans- 
ducer. There  was  a total  of  78  thermocouples,  eight  heat  flow  meters, 
two  solar  pyronometers , plus  transducers  for  wind  speed,  wind  direction 
and  rainfall.  A special  transducer  for  ultraviolet  radiation  data  was 
used  for  a limited  time  period  for  the  Building  Materials  Section  of 
NBS . This  Section  operates  an  outdoor  weatherabi lity  study  station  in 
the  vicinity  of  the  earth  temperature  station. 
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7.1  Temperatures 


The  majority  of  the  78  thermocouples  were  for  soil  temperature 
measurements.  Thermocouples  were  located  at  depths  of  30',  20',  10', 

6',  4',  2',  1',  12",  4",  1",  and  0".  The  12"  depth  temperatures  were 
obtained  using  two  thermocouples. 

For  Patch  #3  the  temperatures  at  depths  of  30’,  20',  10',  6',  4', 

2',  and  1'  were  repeated  using  channels  37  through  43  in  addition  to  the 
regular  thermocouples  which  occupied  channels  25  through  36.  The  seven 
additional  thermocouples  installed  in  patch  #3  were  unprotected  in  the 
earth  along  the  outside  of  the  probe  to  study  the  effect  of  ground  loop 
electric  leakage  (common  mode  noise). 

Thermocouple  channels  68  through  75  were  designated  for  temperature 
measurements  in  and  around  the  underground  instrument  bunker.  Because 
the  bunker  houses  electrical  instruments  which  generate  heat,  and  since 
the  thermal  environment  of  the  instrument  bunker  itself  is  important 
the  thermocouples  were  placed  to  measure  the  inside  surface  temperatures, 
ground  temperature  surrounding  the  bunker,  and  air  temperature  above 
it  (Figure  8). 
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Channels  76  and  77  were  assigned  for  measurement  of  dry-bulb  and 
wet-bulb  temperature  of  air  unaffected  by  radiation  from  the  sky  and  the 
ground.  As  shown  in  Figure  8 these  two  thermocouples  were  installed  in 
a ventilated  birdhouse- like  enclosure  4'  above  the  ground  (standard 
Weather  Bureau  Instrumentation),  where  the  outdoor  air  is  force-ventila- 
ted through  the  enclosure  using  a small  fan.  The  wet-bulb  thermocouple 
was  attached  to  a metal  washer  1/8"  in  diameter  which  was  inserted  into 
a wet  wick.  Capillary  action  kept  the  wick  wet  continuously  because 
the  other  end  of  the  wick  was  in  a plastic  bottle  containing  distilled 
water.  This  scheme  was  recommended  in  reference  [7].  The  thermocouple 
connected  to  Channel  78  was  for  measurement  of  the  temperature  at  a 12' 
height  where  wind  transducers  were  installed  (Figure  8). 

7.2  Solar  Radiation 

Channels  79  through  82  were  assigned  for  pyronometers  (180°  total 
solar  radiation)  to  measure  direct  and  diffuse  radiation.  Only  channels 
80  and  81  were  used  and  these  were  connected  to  shaded  and  unshaded 
pyronometers  (Epply  type),  respectively.  The  calibration  of  the 
pyronometers  for  the  shaded  application  was  done  at  the  Epply  laboratory 
whereas  the  pyronometer  for  the  unshaded  application  was  calibrated  at 
the  ESSA  laboratory.  Solar  radiation  is  very  important  in  the  evalua- 
tion of  the  ground  surface  effect  upon  the  earth  temperature  beneath  it. 


15 


Since  the  temperature  measuring  site  had  no  natural  shade,  the  total 
hemispherical  radiation  measured  by  the  pyronometer  is  considered  most 
meaningful.  Shaded  pyronometer  measurements  were  collected  simultaneously 
to  obtain  data  on  diffuse  sky  radiation  and  to  check  theory  relating 
direct  radiation  and  diffuse  sky  radiation.  Direct  radiation  is  the 
difference  between  the  total  hemispherical  and  the  diffuse  sky  radia- 
tion . 

Since  the  total  hemispherical  radiation  fluctuates  very  rapidly 
during  a cloudy  day,  particularly  when  the  clouds  are  moving  rapidly 
between  the  sun  and  the  sensor,  a once-an-hour  observation  of  the  in- 
stantaneous reading  of  the  pyronometer  is  not  as  meaningful  as  an 
average  reading  that  is  based  upon  an  integrated  measurement  over  a 
one  hour  period.  Such  a long-time  integrating  device  was  not  available 
at  the  start  of  this  experiment  but  has  just  recently  become  available. 

It  is  recommended  that  for  future  data  logging  operations  of  this  nature 
a long-time  integrator  be  used  to  obtain  what  is  essentially  fluctuating 
data.  This  would  apply  to  the  variables  of  solar  radiation,  wind  ve- 
locity, heat  flux  and  possibly  others. 
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7.3  The  Heat  Flow  Data 


Channels  85  through  92  were  connected  to  heat  flow  transducers 
obtained  from  the  National  Instrument  Laboratory,  Inc.  type  no.  HF3 . 
These  transducers  are  polyvinyl  discs  0.16"  thick  and  4.3"  in  diameter 
and  were  considered  suitable  for  the  underground  installation.  The 
manufacturer's  calibration  constants  were  used  for  the  application  of 
these  heat  flow  meters. 

Except  for  patch  #3,  all  the  heat  flow  meters  were  installed  at 
depth  of  4".  Heat  flow  was  measured  at  18",  12",  4"  and  2"  depths  for 
patch  #3  to  study  the  effect  of  earth  cover  as  combined  with  other 
parameters  upon  the  ground  heat  transfer. 

7.4  Wind  Data 

A wind  direction  and  velocity  transducer,  type  C,  which  was  made 
by  the  Belfort  Instrument  Company,  was  assigned  to  channels  83  and  84 
respectively.  The  wind  velocity  transducer  is  of  the  dynamometer  type 
and  its  DC  voltage  outputs  were  calibrated  against  known  air  velocities 
in  the  wind  tunnel  of  the  National  Bureau  of  Standards.  The  calibration 
curve  is  shown  in  Figure  9.  This  transducer  was  installed  atop  a 12' 
tower  as  shown  in  Figure  8. 

The  wind  direction  transducer  is  of  the  vane-type.  The  particular 
transducer  employed  for  this  study  was  so  designed  to  make  eight  switch 
contacts  at  the  base  of  the  van  shaft  corresponding  to  the  eight  wind 
directions;  N,  NE,  E,  SE,  S,  SW,  W and  NW . 
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A signal  conditioner,  shown  in  Figure  10,  was  designed  and  con- 
structed to  generate  eight  different  DC  voltage  levels  corresponding 
to  the  above  eight  wind  directions.  Table  II  shows  the  calibration 
data  as  voltages  vs.  wind  direction. 

7.5  Rain  Gauge 

Channels  94  and  95  were  assigned  for  collection  of  rainfall  data. 
The  rainfall  transducer  was  made  by  Belfort  Instrument  Company,  Catalog 
No.  5915  (Figure  11).  It  consists  of  a galvanized  bucket  placed  on  a 
weighing  mechanism,  which  converts  the  weight  gain  of  the  bucket  due 
to  accumulated  rain  into  a change  of  electrical  resistance.  Figure  12 
shows  the  calibration  curve  provided  by  the  manufacturer  relating  the 
electrical  resistance  to  inches  of  rainfall.  Before  use,  the  calibra- 
tion curve  was  checked  for  its  accuracy  at  the  points  indicated  by 
solid  triangles.  Since,  as  mentioned  later,  the  data  acquisition 
system  employed  for  this  project  did  not  have  the  capability  of  reading 
the  electrical  resistance  directly,  a signal  conditioner  to  convert 
the  resistance  change  to  a mV  signal,  Figure  13,  was  devised.  In 
Figure  13  the  equation  indicated  was  used  to  obtain  the  number  of 
inches  of  rainfall  from  the  DC  signals  sent  to  channels  94  and  95. 
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7.6  Soil  Moisture  Transducers 


Although  not  connected  to  the  automatic  data  acquisition  system, 
the  soil  moisture  content  was  monitored  concurrently  with  the  collection 
of  other  data  but  readings  were  taken  only  once  a week  because  change 
of  moisture  content  was  extremely  slow.  The  moisture  content  trans- 
ducers (Soil  Test  Inc.,  MC310A)  were  installed  at  the  depths  of  4"  and 
12"  at  all  patches  except  Patch  #4  where  the  soil  moisture  gages  were 
located  at  depths  of  30',  20',  10' , 6',  4',  2',  and  1*.  The  moisture 
sensor  was  made  of  parallel  sheets  of  electrode  mesh  separated  by  a 
thin  glass  fiber  cloth  (Figure  14).  This  thin  material  very  quickly 
establishes  a moisture  equilibrium  with  the  surrounding  soil.  When 
an  AC  voltage  is  applied  to  the  electrode  current  leaks  through  the 
glass  fiber  cloth  in  proportion  to  its  moisture  content.  These  gages 
require  calibration  for  the  specific  type  of  soil  in  which  the  gauges 
are  to  be  used.  Thus  a calibration  curve  of  AC  current  versus  a 
known  soil  moisture  content  must  be  obtained.  NBS  laboratory  tests 
using  the  soil  of  the  temperature  measuring  site  show  reasonably  good 
correlation  between  the  moisture  meter  reading  and  the  soil  moisture 
content . 


19 


There  are  many  other  types  of  soil  moisture  gages.  It  was  found 
that  the  particular  transducer  employed  in  this  project  was  satisfactory 
especially  from  the  standpoint  of  use  with  an  automatic  data  acquisition 
system.  A primary  requirement  for  a soil  moisture  gage  was  that  it 
should  be  fast  in  response  yet  remain  very  stable  over  a long  period 
of  time,  and  also  that  its  signal  output  must  be  compatible  with  the 
automatic  data  acquisition  system. 

8.  Data  Acquisition  and  Processing  System 
8.1  Overall  System 

Figure  16  shows  the  overall  data  flow  diagram  of  the  acquisition 
and  data  processing  system  employed  in  this  project.  The  data 
acquisition  system  was  assembled  by  Keltec  Industries  on  contract  and 
consists  of  a scanner,  a digital  voltmeter,  and  a TWX  coupler.  The 
NBS  specification  for  this  system  is  attached  as  Appendix  1 of  this 
report.  Signals  from  the  several  transducers  were  fed  to  the  scanner 
as  analog  signals  (i.e.,  not  in  the  form  of  digits).  Transducers  were 
automatically  scanned  once  an  hour  using  a Cunningham  cross-bar  scanner, 
309G.  The  analog  signals  were  then  digitized  using  a Vidar,  integrating 
type  digital  voltmeter  Model  No.  520.  The  integrating  digital  voltmeter 
(IDVM)  keeps  its  reading  gate  open  for  a period  of  one  second  and  then 
creates  a digitized  mV  value  corresponding  to  the  1 second  average  of 
the  signal  from  the  channel  being  scanned.  The  digitized  data  thus 
created  is  said  to  be  of  a parallel  form  because  the  numerical  values 
converted  into  binary  digits  are  presented  in  a string  of  parallel  bits. 
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Since  the  parallel  bits  are  in  a BCD*  format  they  are  then  converted 
into  a serial  format  by  the  Keltec  TWX  coupler.  This  is  done  to  allow 
the  data  to  be  transmitted  via  a telephone  line  to  a recorder,  bit  by 
bit . 

The  data  as  a string  of  serial  bits  is  transformed  using  a trans- 
mitting data  set  into  tones  which  modulate  the  base  frequency  of  the 
telephone  line  (a  standard  technique  used  by  the  telephone  company). 

The  frequency  modulated  telephone  tone  signal  travels  via  an  underground 
telephone  line  to  a receiving  data  set  which  demodulates  the  message  to 
the  original  serialized  digital  format.  The  signal  then  goes  to  tele- 
type unit,  33TU,  located  in  Room  B307,  Building  226,  of  NBS . The  tele- 
type unit  prints  out  the  data  and  simultaneously  punches  the  same  data 
on  a reel  of  paper  tape.  The  data  on  punched  paper  tape  can  be  read  onto 
a magnetic  tape  later  using  the  NBS  Digidata  1720  paper  tape  to  magnetic 
tape  device.  This  conversion  from  paper  tape  to  magnetic  tape  procedure 
could  have  been  eliminated,  had  the  NBS  computer  system,  UNIVAC  1108, 
been  able  to  read  the  punch  paper  tape.  Unfortunately  this  capability 
was  not  available.  The  magnetic  tape  thus  produced  is  then  read,  bit 
by  bit,  with  the  use  of  a special  computer  program  called  CHAR(N). 


BCD  (Binary  Coded  Decimal)  will  be  explained  in  detail  later. 

Listings  of  all  of  the  computer  programs  mentioned  hereafter  are 
attached  at  the  end  of  this  report. 
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Another,  and  a major  computer  program  called  SOIL  then  reconstructs 
the  original  DC  voltage  data  and  converts  it  into  meaningful  engineering 
information.  If  desired,  the  data  thus  processed  through  SOIL  may  be 
presented  in  graphic  form.  The  processed  engineering  data  are  also 
stored  on  magnetic  tape  for  further  mathematical  and  statistical  analyses. 
These  analyses  will  be  the  subject  of  a separate  report. 

8.2  Scanner  and  Its  Control 

The  scanning  or  multiplexing  of  the  transducers  described  in  pre- 
vious sections  was  accomplished  using  a Cunningham  cross-bar  type 
scanner.  The  cross-bar  is  a system  of  relays  designed  to  open  and 
close  the  switches  rapidly  in  a desired  sequence.  The  scanner  is 
capable  of  switching  600  single  wire  circuits,  which  can  be  made  into 
three  different  modes,  namely,  the  100-6  pole,  200-3  pole  circuits, 

300-2  pole  circuits  or  600-1  pole  circuits.  The  intermixing  of  these 
four  modes  is  usually  difficult.  Since  the  earth  temperature  study 
project  was  concerned  with  about  100  channels,  the  6-pole  mode  was 
adopted;  but  only  three  wires  per  channel  were  used. 

Figure  3 shows  how  the  three  wire  terminals  were  connected  to 
the  various  transducers. 

When  selecting  a scanner,  the  following  specifications  are 
pertinent : 
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1.  Open  circuit  insulation  between  switching  conductors: 
more  than  10  K Q at  0-507,  RH,  -10  °C  - +40  °C  at 

500  VDC 

2.  Maximum  contact  resistance:  less  than  0.3  ohm 

3.  Breakdown  voltage:  1000  VAC 

4.  Current  carrying  capacity  for  non-inducing  load: 

100  mA  at  50  VDC 

5.  Channel  crosstalk:  less  than  60  db  at  5 mega  hz 
with  100  ohm  circuit  impedance 

The  scanner  is  required  to  possess  very  high  leakage  resistance 
and  low  shunt  capacitance  in  order  to  allow  high  source  impedance  and/or 
high  frequency  input  to  be  switched  without  the  degredation  of  the  trans- 
ducer signals. 

It  is  recommended  that  each  channel  be  individually  guarded  and 
floated*  so  that  the  capacitance  between  any  guard  and  the  measurement 
pair  be  extremely  small  (less  than  1 pf).  The  crosstalk  noise,  men- 
tioned in  the  specification,  is  the  spurious  voltage  reading  superimposed 
on  the  source  signal  circuit.  The  noise  is  the  leakage  across  two 
channels  and  is  usually  expressed  in  db,  based  upon  the  following 
equation 

Crosstalk  noise  = 20  log-^  2nfCR 

where  f,  C and  R are  noise  frequency,  interchannel  capacitance  and 
source  impedance,  respectively. 


None  of  the  parts  of  the  circuit  are  grounded. 
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8.3  Integrating  Digital  Voltmeter 


The  heart  of  the  automatic  data  acquisition  system  is  the  inte- 
grating digital  voltmeter  (IDVM),  which  accepts  the  analog  signals  and 
digitizes  them  for  display  as  well  as  for  further  automatic  processing. 

Although  detailed  explanations  of  the  IDVM  operation  can  be  found 
elsewhere  it  is  important  that  the  environmental  engineer  understand 
some  of  the  unique  features  of  its  operation.  Most  of  the  IDVMs  com- 
bine a voltage-to-frequency  (V-F)  converter  with  an  electronic  counter. 
The  V-F  converter  produces  electric  pulses  at  a rate  proportional  to 
the  analog  signal.  The  counter  then  totals  the  pulses  for  a selected 
period  of  time,  providing  a reading  equal  to  an  integrated  average  of 
the  input  during  the  measurement  period.  This  integrating  capability 
helps  to  decrease  error  due  to  superimposed  AC  noise  when  measuring  a 
DC  signal.  When  the  AC  noise  ripple  becomes  large,  however,  the  peak 
value  of  the  ripple  signal  goes  through  zero  and  requires  the  electronic 
counter  to  count  down  (Figure  17).  Many  inexpensive  integrating  digital 
voltmeters  do  not  have  the  capability  of  counting  down  as  well  as  up. 
When  the  use  of  the  up-and-down  counter  is  not  possible,  the  input  AC 
ripple  noise  should  be  made  sufficiently  small  by  other  means  such  as 
guarding.  This  will  be  described  later.  IDVM  employed  in  this  research 
had  no  up-and-down  counter  and  the  counter  frequently  became  overloaded 
when  AC  noise  was  superimposed  on  the  channel  signal  the  magnitude  of 
the  noise  being  much  larger  than  the  DC  signal  being  measured. 
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The  important  features  to  be  considered  when  selecting  an  IDVM  for 
earth  temperature  measurement  are,  its  ability  to  read  low  level  signals 
of  less  than  5 mV , a high  AC  noise  rejection  capability,  a high  input 
impedance,  accuracy  and  resolution,  and  stability  with  respect  to 
ambient  temperature  and  time. 

In  order  to  understand  the  requirements  for  the  proper  functioning 
of  the  IDVM,  it  is  advantageous  to  know  some  basic  characteristics  of 
voltmeters  to  be  used  for  measuring  low-level  signals,  such  as  the  out- 
put of  thermocouples.  There  are  two  important  kinds  of  voltage  measuring 
devices;  one  is  a bridge  balance  or  potentiometric  type  and  the  other  is 
the  current  drawing  type.  The  potentiometric  type  is  based  upon  the 
null -detect ion  of  a galvanometer  in  a Wheatstone  Bridge  circuit,  as 
shown  in  Figure  18.  When  the  bridge  circuit  is  completely  balanced, 
there  is  no  electric  current  flowing  across  the  transducer  circuit. 

As  will  be  described  later,  most  voltage  measurement  errors  come  from 
current  flowing  across  the  transducer  circuit.  The  potentiometric 
voltmeter,  therefore,  is  regarded  as  the  most  accurate  type  of  voltmeter 
when  a stable  source  DC  signal  is  measured.  Unfortunately,  however, 
the  bridge  balancing  and  null-detection  does  not  lend  itself  to  high 
speed,  automatic  data  acquisition,  although  high  speed  is  not  a require- 
ment for  this  particular  investigation.  Most  of  the  automatic  voltage 
recording  devices,  both  analog  and  digital,  are  current -drawing  types, 
or  they  rely  on  a minute  current  in  the  transducer  circuit  to  drive  the 
meter.  There  are  many  cases  where  the  current  drawn  by  the  voltmeter 
presents  enough  voltage  drop  in  the  transducer  circuit  to  cause  sub- 
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stantial  error.  Figure  19  shows  a schematic  diagram  of  the  non- 
potentiometric  voltmeter  circuit  and  indicates  an  equation  relating  the 
actual  voltage  Vs,  developed  by  the  thermocouple  to  V,  the  voltage  in- 
dicated on  the  voltmeter  when  the  transducer  circuit  impedance  and  the 
voltmeter  impedance  (input  impedance)  are  Rg  and  R^,  respectively. 

In  order  for  V and  to  be  very  close,  the  input  impedance  R^  of  the 
non-potentiometric  voltmeter  should  be  many  times  greater  than  the 

source  or  transducer  circuit  impedance  R (R  /R  0). 

s s n 

For  the  earth  temperature  instrumentation,  the  thermocouple  cir- 
cuits are  usually  small  gauge  (#24)  and  very  long  (600  ft)  resulting 
in  a large  value  of  R^ . The  measured  thermocouple  circuit  resistance 
for  this  project  was  on  the  average- 300  ohms.  The  input  impedance  of 
the  voltmeter  should,  therefore,  be  not  less  than  30  K,  if  1%  accuracy 
is  desired.  If  the  voltmeter  has  a relatively  small  input  impedance 
a presimp lifer  should  be  used  to  increase  the  input  impedance  of  the 
system. 


8.4  Noise  Rejection  of  Voltmeter 

During  the  selection  of  a digital  voltmeter,  environmental  en- 
gineers are  often  confronted  with  unfamiliar  terminology  such  as 
common  mode  and  normal  mode  rejection.  And  these  terms  are  expressed 
in  decibels  (db)!  It  is  important  that  these  be  correctly  understood. 
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The  common  mode  noise  is  generated  by  the  AC  line  pickup,  or  hum, 
resulting  from  having  two  grounds  in  the  system;  one  at  the  transducer 
and  another  at  the  chassis  of  the  voltmeter  (Figure  20b).  The  voltage 
difference,  e between  these  two  grounds  induces  a ground  loop  currenL 
at  the  AC  line  frequency  and  causes  a voltage  to  appear  at  the  voltmeter 
input  which  can  easily  be  greater  than  the  signal  itself,  yielding  a 
completely  erroneous  reading.  The  degree  to  which  this  ground  loop 
current  is  minimized  depends  upon  the  electric  isolation  between  the 
transducer  circuit  and  the  voltmeter  chassis.  Use  of  a heavy  ground 
bus  or  shielding  of  the  transducer  cable  is  usually  found  ineffective 
for  decreasing  the  common  mode  pick  up,  and  sometimes  it  is  even  found 
to  increase  it.  Most  of  the  modern  digital  voltmeters  are,  however, 
of  guarded  construction  and  have  a very  high  common  mode  noise  rejec- 
tion capability. 

Figure  21  illustrates  unguarded  and  guarded  voltmeters  and  the 
manner  in  which  ground  loop  current  flows  through  the  transducer  cir- 
cuit. Included  in  the  same  figure  is  mathematical  formula  for  evalu- 
ating the  common  mode  noise  rejection  capability  as  a function  of  noise 
frequency  f,  capacitance,  C1  between  the  transducer  circuit  and  the 
voltmeter  chassis,  and  the  impendance  R of  the  transducer  circuit. 
According  to  this  equation,  the  noise  rejection  of  160  db  means  that 
the  noise  voltage  induced  in  the  transducer  circuit  due  to  common  mode 
voltage  of  100  V is  on  the  order  of  1 microvolt.  The  noise  rejection 
capability  of  the  voltmeter  decreases  by  20  db  with  the  tenfold 
increase  of  any  one  of  f,  C'  and  R. 
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The  normal  mode  noise,  on  the  other  hand  is  generated  by  electric 
appliances,  radios,  and  TV's.  It  will  be  superimposed  on  the  source 
signal  in  the  transducer  circuit,  as  shown  in  Figure  20A.  A conventional 
way  to  minimize  the  normal  mode  noise  pickup  is  to  shield  the  transducer 
lines.  The  normal  mode  noise  can,  however,  be  rejected  completely  by 
the  integrating  process  when  a whole  number  of  noise  cycles  is  averaged 
over  a given  sampling  or  integrating  period  such  as  illustrated  in 
Figure  22.  If  a portion  of  the  noise  cycle  is  left  out  of  this 
averaging  process  (as  shown  in  Figure  23A),  an  error  will  be  produced. 

A maximum  value  for  this  error  is  indicated  in  the  figure.  The  magnitude 
of  error  thus  generated  will  decrease  as  the  noise  frequency  increases, 
becoming  zero  every  time  the  noise  frequency  becomes  an  exact  multiple 
of  the  sampling  frequency  (Figure  23B).  If  the  sampling  period  is 
0.0166  seconds  (60  Hz),  the  normal  mode  rejection  becomes  infinity  for 
noises  of  60  Hz,  12  0 Hz,  180  Hz,  240  Hz,  etc.  Therefore,  except  at 
this  peak  noise  rejection,  the  normal  mode  noise  rejection  generally 
increases  by  20  db  when  the  frequency  of  the  noise  increases  by  a 
decade,  an  opposite  of  the  common  mode  noise  rejection. 

As  indicated  earlier,  the  common  mode  noise  rejection  capability 
of  a digital  voltmeter  decreases  to  20  db  per  decade  increase  in  noise 
frequency,  resulting  in  the  generation  of  normal  mode  noise  on  the 
transducer  output.  The  normal  mode  noise  rejection  capability  in- 
creases at  the  same  rate.  The  combined  noise  rejection  curve  is 
illustrated  in  Figure  24. 
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8.5  TWX  Coupler 


Before  the  digital  data  produced  in  the  IDVM  are  transmitted  through 
telephone  lines,  the  data  must  be  converted  into  a suitable  form.  The 
TWX  coupler  is  designed  to  accept  data  in  a BCD  (binary  coded  decimal) 
form  and  it  is  presented  to  the  coupler  in  parallel  form,  i.e.,  all 
digits  are  presented  at  the  same  time.  The  TWX  coupler  stores  the 
data,  translates  it  into  ASCII  format  (Figure  25)  and  serializes  it 
for  a bit  by  bit  transmission  through  telephone  lines.  The  coupler 
also  initiates  the  scanning  cycle  for  a preset  time  interval  such  as 
30,  60,  12  0 and  240  minutes. 

8.6  Data  Carrier  Set  (Modem) 

The  digital  data  as  serialized  by  the  TWX  coupler  is  frequency 
modulated  on  a carrier  wave  of  765  Hz  by  a transmission  unit  in  the 
Data  Carrier  Set.  The  signals  converted  into  tones  are  carried  via 
telephone  lines  to  a receiving  unit  of  the  Data  Set  and  demodulated 
to  a digital  format.  The  Data  Set  employed  in  this  project  was 
Stelma  IBTCU/lBRCU  and  had  a data  carrying  speed  of  150  band  (bit 
per  second). 
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Since  each  ASCII  character  is  made  of  8 bits,  the  character  trans- 


mission should  be  approximately  20  per  second.  And  since  the  data  for 
each  channel  consists  of  12  characters,  approximately  2 channel  data 
will  be  transmitted  per  second. 


8.7  Teletype  Unit 


Data  coming  from  the  receiving  unit  of  the  data  set  are  assembled, 
converted  into  decimal  digits,  printed  in  page  format  and  simultaneously 
punched  in  a paper  tape  in  ASCII  code  by  the  teletype  (Model  33)  as 
shown  in  Figure  26. 

The  printout  of  the  teletype  is  also  shown  in  Figure  27.  The 
printing  and  punching  speed  of  the  teletype  unit  is  10  characters  per 
second . 

Each  channel  of  the  100  channels  of  information  printed  or  punched 
out  consist  of  12  characters  in  the  following  sequence 


Channel  identification 


3 characters 


polarity 


1 


It 


data 


6 


II 


decimal  location 


1 


n 


space  or  carriage  return 


II 


For  example,  if  the  channel  information  reads 


060  + 0159677  A,  (A  meaning  a space) 
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the  information  contained  in  channel  60  is 


+ 15967  x 10  ^ v or  1.596  mV 

Since  100  channels  were  scanned  and  recorded  at  every  hour,  the 
total  time  required  to  complete  the  data  recordings  per  cycle  was 

12  x 100/10  = 12  0 second  or  two  minutes 

The  teletype  unit  was  activated  by  a signal  from  the  scanning  system  a 
few  seconds  before  the  data  scanning  started  and  remained  activated 
approximately  one  minute  after  all  the  100  channels  of  data  were  re- 
corded. Timing  control  of  the  teletype  unit  was  accomplished  by  a 
pulse  communication  device  indicated  in  Figure  16. 

The  punch  paper  tape  on  the  teletype  unit  comes  in  a roll  that  can 
record  10,000  sets  of  12  character  data.  Since  each  scanning  cycle 
produces  100  sets  of  12  character  data  at  the  rate  of  one  cycle  per 
hour  one  roll  of  paper  punch  tape  should  be  able  to  record  100  hours 
or  approximately  four  days  of  data. 
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8.8  Processing  of  Paper  Punch  Tape 


The  punch  paper  tape  from  the  teletype  is  fed  to  a paper-tape-to- 
magnetic-tape-converter  (Digidata  1720).  The  tape  converter  reads  the 
paper  tape  at  a rate  of  100  characters  per  second  and  writes  on  a mag- 
netic tape  at  200  characters  per  inch  in  even  parvity.  The  format  of 
the  data  frame  on  the  magnetic  tape  is  very  similar  to  that  of  the 
punch  paper  tape . 

Standard  magnetic  tape  is  usually  2400  feet  long  and  should  be 
able  to  contain  5,760,000  characters  of  character  data. 

Since  one  set  of  data  consists  of  twelve  characters  and  one 
scanning  cycle  contains  100  sets  of  data,  one  magnetic  tape  should  be 
able  to  store  approximately  4800  hours  or  1/2  year  of  100  channel  in- 
formation taken  hourly. 

The  time  required  to  produce  a magnetic  tape  is  16  hours  if 
operated  continuously.  As  mentioned  before,  one  reel  of  paper  punch 
tape  contains  four  days  of  information  and  requires  approximately  20 
minutes  to  convert  into  a magnetic  tape  of  approximately  50  ft.  The 
data  structure  in  the  magnetic  tape  is  170  logical  words  of  36  bits 
length,  made  into  a block.  A typical  block  listing  is  shown  in 
Figure  28. 
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8.9  Data  Processing  on  Computer 


The  magnetic  tape  produced  on  the  Digidata  1720  is  fed  into  a 
computer  (UNIVAC  1108  at  the  National  Bureau  of  Standards)  at  a rate 
of  150  in/sec  or  2500  data/sec.  A computer  program  called  CHAR(N),  the 
Fortran  listing  of  which  is  attached  in  this  report,  reads  the  tape 
block  by  block  and  decodes  the  data,  twelve  binary  characters  at  a 
time,  to  form  BCD  characters.  The  decoding  scheme  of  the  magnetic 
tape  of  CIIAR(N)  is  illustrated  in  Figure  2 9. 

The  data  decoded  by  CHAR(N)  is  processed  through  a computer  program 
called  SOIL  to  obtain  engineering  information.  A Fortran  listing  of 
SOIL  is  attached  to  this  report  and  Figure  30  shows  a sample  listing 
of  one  cycle  of  information  converted  into  engineering  data.  This 
listing  shows  that  channels  26,  41,  and  59  are  not  operating  properly. 
The  columns  headed  by  CHAN,  MV,  DATA  and  REMARK  of  Figure  30  represent 
the  channel  number,  mV  data  acquired  by  the  data  logger,  engineering 
data  converted  from  the  mV  data,  and  appropriate  remarks  respectively. 
The  remarks  indicated  are 

T = temperature  in  degree  F 

2 

SOL  = solar  radiation  in  Btu/hr,  ft 

WDD  = wind  direction 

WDV  = wind  velocity  in  MPH 

2 

HF  = heat  flow  in  Btu/hr,  ft 

RAIN  = rain  data  in  inches  of  water 

REF  = reference  temperatures 

= ultra-violet  radiation 
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UVR 


If  desired,  the  data  processed  through  the  computer  program  SOIL  can  be 
presented  in  graphical  forms  such  as  illustrated  in  Figures  31  and  32. 

9.  Error  Analysis 

In  addition  to  the  stem  conduction  error  of  the  probe  containing 
the  thermocouples,  and  the  noise  errors  discussed  previously,  the 
following  errors  should  also  be  taken  into  consideration  when  using 
thermocouples  for  earth  temperature  measurement. 

9.1  Thermoelectric  Error 

This  is  an  emf  caused  by  having  junctions  of  two  dissimilar  metals 
under  an  appreciable  temperature  gradient  within  the  thermocouple  cir- 
cuit. When  the  junction  of  two  dissimilar  metals  is  unavoidable,  an 
isothermal  environment  should  be  created  around  the  junction. 

9.2  Galvanic  Error 

The  emf  created  by  the  electrolytic  or  galvanic  effect  is  no  small 
problem  for  ground  temperature  measurement.  Thermocouple  circuits 
should  be  kept  dry  and  clean  to  minimize  this  problem,  which  means  the 
complete  protection  of  the  transducer  circuits  from  the  moist  soil. 
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9.3  Electrostatic  and  Triboelectric  Effect 


There  are  low  frequency  disturbances  generated  by  electrostatic 
and  triboelectric  (charge  of  electricity  due  to  friction)  effect.  When 
the  user's  body  potential,  is  quite  different  from  the  ground  potential, 
the  voltmeter  indicates  erroneous  signals  with  the  approach  of  his  hand. 

9.4  Stability  of  the  Voltmeter 

The  voltmeter  accuracy  is  usually  specified  in  +%  of  the  full 
scale  within  a limited  temperature  range  (0  - 50  °C  is  very  common), 
and  for  a certain  time  period.  The  drifting  of  the  voltmeter  reading 
due  to  the  time  and  the  ambient  temperature  change  should  be  carefully 
monitored.  A good  way  to  check  the  voltmeter  stability  is  to  monitor 
its  reference  voltage,  which  is  usually  in  a temperature  controlled 
oven  inside  the  voltmeter.  Channel  000  of  this  system  was  connected 
to  the  reference  source  of  one  volt.  The  reading  was  1.00028  volt  in 
May,  1968  and  has  become  1.00097  in  July,  1969,  which  is  approximately 
+0.005%  drift  per  month  and  thus  considered  extremely  stable.  The 
daily  fluctuation  of  the  reference  source  was  practically  unnoticeable . 
Stability  and  accuracy  of  the  reference  thermocouple  junction  are  also 
important  from  the  standpoint  of  converting  the  thermocouple  emf  to 
actual  degrees  Fahrenheit.  The  reference  temperature  junction  used  in 
this  research  for  channels  1 through  67  was  a 150  °F  source  manufactured 
by  PACE  Company.  The  manufacturer  claimed  that  the  temperature  stability 
was  +1/10  °F  and  junction  accuracy  was  +1/5  °F.  Similar  stability  and 
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accuracy  were  claimed  by  the  Joseph  and  Kays  Company  for  their  32  °F 
reference  junction  used  for  channels  68,  75,  76,  77  and  78. 

Since  these  reference  junctions  were  shared  by  many  thermocouple 
circuits,  it  is  possible  that  the  common  mode  errors  introduced  into 
one  circuit  would  have  influenced  the  others. 

In  order  to  evaluate  the  magnitude  of  the  common  mode  error 
existing  in  the  transducer  system  a comparative  analysis  was  made  of 
the  readings  for  most  of  the  transducers  obtained  by  the  integrated 
digital  voltmeter  and  those  obtained  by  a potent iometric  analog  volt- 
meter (Rubicon  Honeywell  274g).  As  mentioned  earlier,  the  potentiometric 
voltmeter  should  be  devoid  of  all  the  errors  due  to  the  common  and  nor- 
mal mode  noises.  Table  3 shows  the  comparative  readings  of  these  two 
voltmeters  with  the  source  impedance  measured  separately  by  an 
Ohmmeter  (Simpson).  Because  of  the  time  elapsed  between  the  two  sets 
of  measurements  and  because  the  ground  surface  condition  was  rapidly 
changing,  comparisons  of  shallow  depth  thermocouple  data  and  particularly 
the  heat  flow  meters  readings  are  not  meaningful.  The  deep  thermocouples, 
however,  show  an  average  deviation  of  the  IDVM  reading  from  the 
potentiometer  reading  of  approximately  7 pV. 

Table  IV  summarizes  all  of  the  possible  errors  discussed  pre- 
viously. A total  temperature  error  corresponding  to  52p.V  is  2.5  °F. 
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10.  Summary 


Instrumentation  used  in  an  earth  temperature  project  at  the 
National  Bureau  of  Standards  is  described  in  terms  of  its  application 
in  an  automated  data  acquisition  and  processing  system.  The  major 
components  of  the  data  acquisition  system  are  presented  with  particular 
emphasis  on  the  principles  and  use  of  a digital  voltmeter. 

The  several  sources  of  error  involving  thermocouple  circuits 
were  analyzed.  It  was  concluded  that  the  precision  of  temperature 
measurement  probable  in  this  type  of  system  is  about  +2  °F  consistent 
with  previous  experience. 

This  report  is  primarily  addressed  to  the  instrumentation  in  the 
system  but  sample  results  of  earth  temperatures  and  heat  fluxes  are 
included  to  illustrate  the  versatility  of  the  automated  digital  data 
acquisition  system. 

The  data  was  collected,  processed  and  stored  on  magnetic  tapes. 

An  analysis  of  the  data  and  results  is  planned  for  later  publication. 
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Table  I 


Sensor 

Reference  Volt 
Thermocouple 

ii 

99 

91 

91 

99 

II 

81 

If 

II 

88 

II 

II 

II 

98 

II 

II 
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It 
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ft 
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II 

II 
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91 
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II 
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II 

91 

18 
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II 
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19 

81 

II 

II 
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II 
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II 

4 
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II 

It 

81 

II 
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18 

II 

It 

II 

II 
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II 
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II 
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II 

II 
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2 
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II 
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II 
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II 
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II 

II 
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II 
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II 

II 

91 
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II 
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ii 
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2 

II 

II 

It 

ii 

II 

1 

II 
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II 

ii 
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23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 
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Table  1 - Continued 


Sensor 

Thermocouple 

II 

II 

II 

II 

It 

II 

It 

ft 

II 

II 

It 

II 

II 

If 

It 

II 

II 

If 

II 

I f 


Location 


12  in.  under  Patch  No.  2 


4 

2 

1 

0 


it  it 


it  it 


It  IV 


ft  ft 


30  ft  under  Patch  No.  3 

2 Q II  It  II  II  II 

II  II  II  II  II 

£ II  II  II  II  If 

^ II  II  II  II  II 

2 it  ii  ii  ii  ii 

^ ii  ii  ii  ii  ii 

12  in.  under  Patch  No.  3 


4 

2 

1 

0 


II  II 


ft  II 


II  II 


II  II 


30  ft  under  Patch  No.  3 


20 

10 

6 


II  II  II 


II  II  II 


II  II  It 
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42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


Table  1 - Continued 


Sensor  Location 


Thermocouple 

4 

ft 

under 

Patch 

No  • 

3 

II 

2 

It 

IB 

II 

it 

11 

II 

1 

II 

II 

II 

it 

11 

II 

30 

ft 
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No. 

/ 

q- 

II 

20 

II 

ti 

II 

n 

it 

II 

10 

II 

ti 

II 

it 

n 

II 

6 

II 

it 

II 
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n 

II 

4 

II 

ii 

II 

it 

it 

II 

2 

II 

ii 

II 
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ti 

II 

1 

II 

ii 
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it 

ii 

II 
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II 

30 

ft 
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No. 

5 

It 

20 

ii 

II 

ii 

it 

II 

II 

10 

n 

II 

n 

it 

II 

II 
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It 

it 

ti 

II 

II 
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it 

II 

ii 

n 

II 
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Table  1 - Continued 


Channel  No. 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 


Sensor 

Thermocouple 

ii 

it 

f? 

If 

II 

II 

II 

II 

II 

If 

II 

II 


74 


75 


ti 


Location 

2 ft  under  Patch  No.  5 

^ II  II  II  II  II 

12  in • under  Patch  No,  5 

ii  ii  ii  ii 

ii  ti  ii  ii 

it  n ii  it 

II  II  II  M 

Zone  Box  within  the 
Bunker 

Inside  Surface  of  the 
Bunker  Wall 

Inside  Surface  of  the 
Bunker  Ceiling 

Inside  Surface  of  the 
Bunker  Floor 

Outside  Surface  of  the 
Bunker  Wall 

Approximately  2 ft  under 
the  Ground  Surface  and 
4 ft  away  from  the 
Bunker 

Ground  Surface  right 
above  the  Bunker 

2 ft  above  the  Ground 
over  the  Bunker 


4 " 
2 " 
1 " 
0 " 
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76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 


Table  1 - Continued 


Sensor  Location 

Thermocouple  Wet-bulb  Temperature  of 

Outside  Air 

" Dry-bulb  Temperature  of 

Outside  Air 

" Air  Temperature  12  ft 

above  the  Ground 
Surface 

19 

Solar  Pyronometer 

ii 

91 


Wind  Direction 
Windspeed 


Heat flow  Meter 

4 

in. 

under 

Patch 

No.  1 

18 

4 

in. 

under 

Patch 

No.  2 

II 

18 

in. 

under 

Patch 

No.  3 

II 

12 

in . 

under 

Patch 

No.  3 

II 

4 

in. 

under 

Patch 

No.  3 

18 

2 

in. 

under 

Patch 

No.  3 

It 

4 

in. 

under 

Patch 

No.  4 

II 

4 

in. 

under 

Patch 

No.  5 

Thermocouple 

Replacement 

26 

for  Channel 

Rain  Gage 

Above 

the  Bunker 

Rain  Gage 

Above 

the  Bunker 

46 


Table  1 - Continued 


Channel  No. 


Sensor 


96 

Thermocouple 

97 

Ultra-violet 

Radiation 

98 

99 

Locat ion 

Chassis  of  the  Data 
Logger 

Above  the  Bunker 


47 


Table  II  Wind  Direction  Conditioner  Output 


Direction 

Output  in  mV 

North 

.31 

Northwest 

.49 

West 

.73 

Southwest 

1.03 

South 

1.18 

East 

5.54 

Northeast 

3.51 

48 


2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 


Table  III  Potentiometric  Check  of  the 


Keltec  Data  Acquisition  System 


Patch 

Depth 

in. 

Resis 

tance 

ohms 

1 

360 

190 

1 

240 

190 

1 

12  0 

180 

1 

72 

180 

1 

48 

180 

1 

24 

180 

1 

12 

190 

1 

12 

190 

1 

4 

190 

1 

2 

190 

1 

1 

190 

1 

0 

190 

2 

360 

170 

2 

240 

160 

2 

120 

155 

2 

72 

155 

2 

48 

155 

2 

24 

155 

2 

12 

160 

Rubicon 

mV 

IDVM 

mV 

Rubicon 

IDVM 
micro  V 

-2  .232 

-2.248 

16 

-2 . 166 

-2 .176 

10 

-1.870 

-1.892 

22 

-1.692 

-1.696 

4 

-1.594 

-1.598 

4 

-1.527 

-1.534 

7 

-1.518 

-1.517 

-1 

-1.511 

-1.514 

3 

-1.442 

-1.425 

-17 

-1.195 

-1.235 

40 

-1.229 

-1.265 

36 

-1.052 

-0.982 

70 

-2.255 

-2  .261 

6 

-2.179 

-2 . 182 

3 

-1.900 

-1.910 

10 

-1.871 

-1.877 

6 

-1.867 

-1.874 

7 

-1.887 

-1.894 

7 

-1.901 

-1.901 

0 
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Channel 

Patch 

Depth 

in. 

20 

2 

12 

21 

2 

4 

22 

2 

2 

23 

2 

1 

24 

2 

0 

25 

3 

360 

26 

3 

240 

27 

3 

120 

28 

3 

72 

29 

3 

48 

30 

3 

24 

31 

3 

12 

32 

3 

12 

33 

3 

4 

34 

3 

2 

35 

3 

1 

36 

3 

0 

37 

3' 

360 

38 

3 * 

240 

39 

3' 

120 

40 

3 9 

72 

41 

3 1 

48 

Rubicon 

IDVM 

Rubicon 

IDVM 

mV 

mV 

micro  V 

-1.946 

-1.954 

8 

-1.841 

-1.836 

-5 

-1.790 

-1.791 

1 

-1.758 

-1.763 

5 

-1.658 

-1.631 

-2  7 

-2.270 

-2.281 

11 

— 

— 

-- 

-1.810 

-1.818 

8 

-1.920 

-1.928 

8 

-1.844 

-1.851 

7 

-1.884 

-1.890 

6 

-1.819 

-1.826 

7 

-2.028 

-2.040 

12 

-1.653 

-1.646 

-7 

-1.593 

-1.603 

10 

-1.544 

-1.551 

7 

-1.492 

-1.179 

-2 13 

-2  .246 

-2  .239 

-7 

-2.215 

-2.219 

4 

-2.028 

-2.041 

13 

-1.926 

-1.935 

9 
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165 
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240 

220 

230 

220 

230 

220 

230 
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260 

240 

240 

240 
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42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 


itch 

Depth 

in. 

Resis- 

tance 

ohms 

3 ' 

24 

240 

3' 

12 

240 

4' 

360 

210 

4' 

240 

210 
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120 

200 
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72 

2 00 

4 

48 

2 00 

4 

24 

200 

4 

12 

200 

4 

12 

200 

4 

4 

2 05 

4 

2 

2 00 

4 

1 

2 00 

4 

0 

2 00 

5 
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300 

5 

240 

300 

5 

120 

300 

5 

72 

300 

5 

48 

300 

5 

24 

2 80 
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300 

5 
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300 

5 

4 

2 90 

Rubicon 


Rubicon 

IDVM 

IDVM 

mV 

mV 

m 1 1 n i 

- 1 . 8-»  3 

- 1 . 8 V) 

13 

-1.818 

-1.823 

5 

-2.262 

-2  .269 

7 

-2.259 

-2.266 

7 

-2 . 12  9 

-2.137 

8 

-2.023 

-2.030 

7 

-1.962 

-1.70 

8 

-1.925 

-1.931 

6 

-1.913 

-1.919 

6 

-1.901 

-1.904 

3 

-1.823 

-1.834 

11 

-1.799 

-1.801 

2 

-1.767 

-1.768 

1 

-1.644 

-1.486 

-178 

-2  .222 

-2  .230 

8 

-2.218 

-2  .233 

5 

-1.892 

-1.898 

6 

-1.919 

-1.928 

9 

-1.899 

-1.902 

3 

-1.888 

-1.898 

10 

-2.072 

-2.078 

6 

-1.837 

-1.844 

7 
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Channel  Patch  Depth 

in. 

Resis- 

tance 

ohms 

Rubicon 

mV 

IDVM 

mV 

Rubicon 
IDVM 
micro  V 

65 

5 2 

2 90 

-1.822 

-1.833 

11 

66 

5 1 

300 

-1.840 

-1.858 

18 

67 

5 0 

300 

-1.451 

- 1 . 14 1 

-10 

68 

Zone  Box 

+1.170 

1.168 

2 

69 

Bunker  Wall 

70 

" Ceiling 

-0.085 

-0.089 

4 

71 

" Floor 

-0.072 

-0.076 

4 

72 

" Soil 

-0.109 

-0.112 

3 

73 

Soil 

-0.256 

-0.260 

4 

74 

Earth  Surface 

+0.073 

+0.077 

4 

75 

2 ' above  Ground 

+1.166 

+1.221 

-55 

76 

DB 

+1.120 

+1.160 

-40 

77 

WB 

+1.113 

+1. 184 

-71 

78 

12  1 above  Ground 

+1.151 

+1.215 

-64 

79 

80 

Solar 

+1.368 

+2 .353 

-985 

81 

82 

83 

Wind  air 

84 

Wind  Speed 

+0.038 

+0.02  9 

9 

85 

Heat  Flux 

-6.800 

-6.197 

-603 

86 

Heat  Flux 

-3.885 

-3.701 

-184 

87 

Heat  Flux 

+0.764 
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Channel  Patch 

Depth 

in. 

Resis- 

tance 

ohms 

Rubicon 

mV 

IDVM 

mV 

Rubicon 

IDVM 

micro  V 

88 

Heat  Flux 

+0.377 

0.537 

-160 

89 

Heat  Flux 

90 

Heat  Flux 

+8.334 

10.324 

-1990 

91 

Heat  Flux 

-3.005 

2.909 

-96 

92 

Heat  Flux 

-3.4  72 

-3.340 

-132 

93 

-2.233 

-2  .22  7 

-6 

94 

Rain  Gage 

-0.063 

-0.062 

-1 

95 

Rain  Gage 

-0.002 

96 

Bunker  Air 

+2 .756 

2.751 

5 

97 

98 

99 
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Table  IV  Summary  of  Temperature  Measurement  Errors 


Thermoelectric  error 

5 microvolts 

Galvanic  error 

2 microvolts 

Triboe lectric  error 

1 microvolt 

Drift  or  stability  error 

10  microvolts 

A 

Normal  mode  error 

> 

Common  mode  error 

7 microvolts 

Stem  conduction  error 

2 microvolts 

Reference  junction  error  (0.5  °F) 

10  microvolts 

Thermocouple  wire  calibration  (0.75  °F) 

5 microvolts 

Total 


42  microvolts 
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The  National  Bureau  of  Standards  is  located  near  Gaithersburg,  Maryland,  just  cff 
Interstate  Route  70S,  20  air  miles  from  the  zero  milestone  in  Washington.  The  exits 
from  the  highway  are  clearly  marked  "National  Bureau  of  Standards." 

Visitors  should  go  to  the  receptionist  in  the  main  lobby  of  the  Administration  Build- 
ing for  directions  or  for  shuttle  service  to  other  buildings  on  the  site. 


rmcNDSHe  I 
NTCNNATKWAL 
AIAPOAT^  f 


DEPARTMENT  OF  COMMERCE 
15th  STREET  AND 
CONSTITUTION  AVENUE,  N.W. 


SRADJTOCk  roa 

CLEARINGHOUSE 


DRIVING  INSTRUCTIONS: 

From  Downtown  Washington  - 

Connecticut  Avenue  to  the  Beltway  and  Route  70S 
North;  or 

River  Road  (Route  190)  to  the  Beltway,  or 
George  Washington  Memorial  Parkway  to  the  Belt- 
way 

Beltway  to  Route  270  and  Route  70S  North 

From  Dulles  Airport  — 

Dulles  Access  Highway  to  Beltway. 

Beltway  to  Route  270  and  Route  70S  North 

From  Notional  Airport  — 

George  Washington  Memorial  Parkway  to  the  Belt- 
wav 

Beltway  to  Route  270  and  Route  70S  North 
From  Friendship  Airport  — 

Washington  - Baltimore  Parkway  to  Beltway 
Beltway  to  Route  70S  North 


OTHER  TRANSPORTATION: 

There  is  no  satisfactory  public  transportation  from  the 
Washington  area  to  the  Bureau  at  Gaithersburg.  How- 
ever the  Bureau  provides  a shuttle  service  for  official 
visitors  and  staff  leaving  the  Department  of  Commerce 
(15th  Street  entrance  nearest  Constitution  Avenue) 
every  hour  8:10  a.m.  to  4:10  p.m.  Travel  time  is  approx- 
imately 40  minutes. 

Return  service  on  the  same  schedule  leaves  from  the 
front  of  the  Administration  Building. 

Please  make  arrangements  in  advance  for  a seat  on  the 

shuttle.  Call  921-2615. 


Note:  The  George  Washington  Memorial  Parkway  is  on 

the  Virginia  side  of  the  Potomac. 
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Figure  2 Detailed  Topographic  Map  of  the  Earth  Temperature  Site 
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Installation  and  Other  Related  Figures 
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Patch  *1  2 
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Figure  3B  Connection  of  Various  Transducers  to  the  Scanner  Terminals 


i 


60 


Figure  3C  Earth  Temperature  Probe  and  Transducer  Conduit 
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Figure  3D  Assignment  of  Heat  Flow  Meters  and  Moisture  Gauges 


Figure  4 Laboratory  Model  of  a Thermocouple  Earth  Temperature  Probe 
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115  V 

250  WATTS 
INFRARED  LAMPS 


Error  Analysis 
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Stainless  Steel  Probe  Test 
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Figure  6 Comparison  of  Sand  Temperature  Measured  by  a Vertical  Probe  and  by  Thermocouples 
Horizontally  Installed  in  the  Sand  at  the  Same  Depth 
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Figure  6A  Comparison  of  Sand  Temperature  Measured  by  a Vertical  Probe  and  by  Thermocouples 
Horizontally  Installed  in  the  Sand  at  the  Same  Depth 


Interior  Dimension  of  Instrument  Bunker 
8'  L,  4'W,  And  7'  H 


FACILITIES  IN  THE  UNDERGROUND  INSTRUMENT  BUNKER 


Figure  7 Underground  Instrument  Bunker  and  Its  Facilities 
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X Thermocouple  Locations  in  and  around  Instrument  Bunker 


Figure  8 Thermocouple  Installation  In  and  Around  the  Bunker 


67 


150 

140 

130 

120 

110 

100 

90 

80 

70 

60 

50 

40 

30 

20 

10 


RPM 

100  200  300  400  500  600  700  800  900  1000 


m V 


Calibration  Curve  for  the  Wind  Velocity  Gauge 
6 8 


Figure  9 


WIND  DIRECTION  CONDITIONER 


Figure  10  Wind  Direction  Signal  Conditioner  Circuit 
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Figure  11  Rain  Gauge 
70 


RAIN  GAUGE  CALIBRATION 


Oo>eoi^toin<frrocvi“-o 


ZTX  ‘30NVlSIS3d 


RAINFALL,  INCHES 


94  95 


R : A Linear  Function  of  RAIN  FALL  (Inches) 


RAIN  GAUGE  SIGNAL  CONDITIONER 


Figure  13  Rain  Gauge  Signal  Conditioner 
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TYPICAL  CONSTRUCTION  DETAILS 


HALF-CASE 
OF  SOIL  UNIT 

(Two  required  per  unit) 


ELECTRODE  SANDWICH 
WRAPPING  DIAGRAM 
(End  view,  schematic) 


ASSEMBLY  VIEW 

(Electrode  Sandwich  cut  away. 
Top  half  case  removed) 


THERMISTOR 


ASSEMBLED 
SOIL  UNIT 


(half-case  flanges 

SPOT  WELDED  TOGETHER) 


Figure  14  Soil  Moisture  Gauge 
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Elapsed  Days 


ANLOG  SIGNALS  IOO  Channels 


Telephone  | Digital 

| ( Frequency  Modulated) 

J Serial 


OUTPUT 

PRINT  OUT 
PAPER  PUNCH 
TAPE 


KELTEC  CP-2330  SYSTEM 


Figure  16  Overall  Flow  Diagram  for  the  Data  Acquisition  System 
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Continuation  of  Figure  16 
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Voltage 
Indicated  = 0 


Erroneous 

Voltage 

Indicated 


Effect  of  up-down  counter  on  low-level  D. C.  signal 


Figure  17  AC  Noise  Crossing  of  Zero  Line  on  Low  Level  DC  Signal 
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Figure  29  Decoding  Scheme  of  the  Magnetic  Tape  Data 
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Figure  31  4"-Depth  Temperature  Under  Asphalt  Cover  (Patch  #1) 
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Figure  32  4"-Depth  Heat  Flow  Meter  Under  Asphalt  Cover  (Patch  #1) 


Appendix  1 


Specifications  for 

Automatic  Data  Acquisition  System  for  Earth  Temperature  Measurement 
(Contractual  Research  Under  Office  of  Civil  Defense, 

Department  of  the  Army,  Contract  OCD-OS-62-44 , Unit  1211A) 

This  research  contract  calls  for  hourly  measurement  of  undisturbed 
earth  temperatures  at  selected  depths  for  five  sites  during  a period  of 
two  years,  starting  from  late  1967.  Also  included  in  the  study  are  the 
hourly  measurements  of  air  dry-and  wet-bulb  temperatures  at  four  feet 
above  the  ground  surface,  and  solar  incident  radiation  and  heat  flux 
at  the  ground  surface.  Temperature  measurements  will  be  made  by  thermo- 
couples (cold- junction  compensation  supplied  by  NBS)  at  ambient  condi- 
tions . 

It  is  contemplated  at  this  moment  that  the  data  acquisition  system 
will  be  housed  in  an  underground  bunker  and  the  data  acquired  will  be 
transmitted  through  telephone  cable  to  a teletypewriter  in  Room  B307, 
Building  226  of  the  National  Bureau  of  Standards.  The  telephone  cable 
distance  is  approximately  one  mile. 

1.  The  system  shall  consist  of  a low-level  scanner  for 

millivolt  signals  from  copper  oonstantan  thermocouples, 
integrating  digital  voltmeter  (or  an  equivalent  digital 
voltmeter  capable  of  acquiring  statistical  average  of 
one  second  sampling  period),  a coupler  to  the  tele- 
typewriter, and  timer.  The  timer  is  to  control  the 
initiation  of  the  scanning  cycle  as  well  as  to  identify 
the  scanning  cycle  by  providing  identification  signal. 
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The  system  shall  be  capable  of  accepting  two  external 
digital  signals  which  will  be  compatible  in  electri- 
cal characteristics  with  the  output  from  the  digital 
voltmeter.  The  external  system  to  produce  two  digital 
signals  will  be  considered  optional  for  this  specifi- 
cation and  the  detail  pertaining  to  this  sub-system 
is  explained  later*. 

2.  A suggested  sequence  of  the  output  will  be  scan 
identification,  data  from  two  external  digital 
sub-systems,  calibration  data,  followed  by  the 
channel  information  consisting  of  channel  identi- 
fication, polarity,  data,  decimal  or  exponent. 

3.  The  paper  tape  shall  be  of  ASC  II  code  and  shall 
have  the  following  word  format  corresponding  to 
each  of  100  channel  data:  channel  identification; 

3 characters,  polarity;  one  character,  data;  6 

characters,  decimal;  one  character,  and  space; 

) 

one  character. 

4.  Each  two  wire  input  to  the  system  shall  be  in- 
dividually guarded  and  interchannel  capacitance 
shall  not  exceed  0.5  pF.  A random  access  mode 

of  operation  shall  be  provided  within  the  scanner. 
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5.  The  integrating  digital  voltmeter  (IDVM)  shall  have 
a sampling  period  of  either  one  second  or  0.166 
second.  If  it  is  a successive  approximation  type, 
the  digital  voltmeter  shall  be  capable  of  acquiring 
the  best  statistical  estimate  of  average  data  for 
one  second  when  operating  under  this  integrating 
period  of  one  second.  The  overall  effective  common 
mode  rejection  of  the  system  shall  be  greater  than 
120  db  and  the  normal  mode  rejection  shall  be  greater 
than  35  db  for  all  the  frequency  above  59  cps. 

6.  Since  the  system  has  to  be  operated  unattended  for 
a period  of  several  weeks  between  the  consecutive 
calibration  checks,  it  is  mandatory  that  it  possess 
a small  drift  characteristic  due  to  time  and  due  to 
ambient  temperature  fluctuations.  The  input  signal 
to  the  scanner  is  low  level,  (a  few  microvolts  to  a 
maximum  of  10  millivolts).  With  respect  to  this 
input  signal  level,  the  required  resolution,  time 
stability  and  temperature  stability  shall  be 
specified  as  follows: 
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System  Resolution  0.0017o  or  better  (0.1 

microvolt  for  10  milli- 
volt signal  level). 

Time  Stability  0.05%  per  week  or  better 

0.17,  per  month  or  better 

Temperature  Stability  + 0.107o  per  deg.  C or 

better  (+  10  microvolts 
for  10  deg.  C change  of 
ambient  at  10  millivolt 
level) . 

7.  The  system  shall  include  a + 1 volt  calibration  source 
in  a temperature  controlled  bath.  The  reference  voltage 
shall  have  drift  characteristics  of  less  than  0.0057,  for 
30  days  and  temperature  coefficient  of  0.0037,/C  for  10 
to  50  C.  This  calibration  voltage  shall  be  provided  to 
the  input  of  the  scanner  and  shall  be  recorded  as  the 
first  channel  of  information  for  each  scanning  cycle. 

8.  The  system  shall  have  an  input  impedance  of  not  less 
than  10  megohms  for  10  MV  signal  range. 

9.  The  time  shall  be  provided  to  allow  manual  selection 
of  time  interval  between  the  scan  initiation.  The 
time  interval  may  be  established  at  30,  60,  120,  or 
240  minutes.  When  the  selected  amount  of  time  has 
passed,  the  timer  issues  an  "initiate"  command  to 
the  scanner,  which  advances  to  the  first  channel  of 
the  100-channel  group.  When  all  the  channels  have 
been  recorded,  the  time  system  remains  idle  until 
the  next  "initiate"  command  is  received.  The  timer 
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system  shall  include  a digital  counter  which  pro- 
duces the  scan  identification  number  (4  digits)  and 
feeds  it  into  the  digital  system. 

10.  The  system  less  the  teletype  unit  shall  be  smaller  than 
36"  high,  36"  wide  and  36"  deep.  The  system  shall  be 
capable  of  normal  operation  for  the  environment  of 

2 0°  + 10  °C  and  70°  + 2 0%  RH. 

11.  The  system  less  the  teletypewriter  shall  be  mounted 
in  a rigid  metal  cabinet  and  the  metal  cabinet  shall 
have  four  eyebolts  at  the  top  and  four  2-1/2"  wheels 
at  the  bottom  for  handling. 

12.  The  vender  shall  inspect  and  test  the  system  at  the 
National  Bureau  of  Standards  and  certify  that  it  is 
in  working  condition  within  the  requirement  of  this 
specification. 

13.  The  vender  shall  warranty  the  entire  system  for  one 
year  from  the  date  of  acceptance  by  the  National 
Bureau  of  Standards.  Service  shall  be  performed 

at  the  NBS  by  a fully  qualified  service  engineer 
to  be  supplied  by  the  vender  within  five  days 
of  notification. 

14.  Vender  shall  provide  comprehensive  instruction 
manuals  including  calibration  and  trouble  shooting 
procedures,  and  complete  parts  lists  * 
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15.  Bidders  are  requested  to  provide  a complete  set  of 
technical  brochures  describing  major  components  of 
the  system,  which  will  be  used  by  NBS  to  evaluate 
the  proposal. 

16.  The  entire  system  shall  be  delivered  to  the  National 
Bureau  of  Standards  (Gaithersburg,  Maryland)  and 
made  operational  within  120  calendar  days  from  the 
date  of  purchase  order. 

*17.  Due  to  the  very  much  fluctuating  nature  of  the  solar 
incident  radiation  and  ground  surface  heat  flux,  it 
is  necessary  to  obtain  hourly  integrated  values.  A 
suggested  method  is  to  convert  millivolt  signal  of 
+ 0-20  linearly  into  frequency  of  0-1000  cpm  and  to 
totalize  the  hourly  frequency  by  an  electronic  bi- 
directional counter  of  at  least  6 digits.  The  counter 
will  count  up  on  positive  signals  and  down  on  negative 
signals.  The  counter  shall  be  automatically  reset  to 
zero,  if  necessary;  at  the  end  of  the  scanning  cycle. 
The  detail  of  this  integrating  system  design  will  be 
left  to  the  vender.  This  long  integrating  system, 
however,  should  be  able  to  maintain  0.1%  linearity 
and  accuracy  of  reading  over  full  range  with  + 10% 
variation  in  line  voltage.  The  timing  of  the  system 
shall  be  such  as  to  allow  the  completion  of  a fixed 
clock  controlled  integration  period  before  the 
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readiation  and  heat  flux  signals  are  recorded.  That 
is  if  the  timer  is  set  on  a 30  minute  cycle  integra- 
tion of  29  minute  values  will  be  provided.  At  60 
minutes,  58  minutes  and  so  on. 
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741 

NQATA 

0000 

I 

014736 

NDAY 

0000 

I 

014742 

NHOUR 

0000 

I 

000377 

NO 

U01 

NT  APE 

0000 

I 

014752 

NTEST 

0000 

I 

014737 

NYR 

0000 

R 

014763 

PV 

D 12 

REMK 

0000 

I 

014726 

SKIP 

0010 

R 

000000 

TCONF 

0000 

R 

000622 

Tf7 

731 

riMEO 

0000 

I 

014730 

TPETST 

0000 

R 

000217 

VI 

0000 

R 

014764 

V? 

7b5 

rtO  V 

0000 

R 

000363 

YR 

c 

c 

L 

c 

C THIS  PROGRAM  READS  DATA  FROM  MAG  TAPE  OR  DATA  CARDS 

C 

C KCARD  = BLANK  SPACE  FOR  CARD  INPUT 

C KCARD  = 1 FOR  TAPE  INPUT 

C SKIP..  NUMBER  OF  BLOCKS  TO  BE  SKIPPED  ON  OUTPUT  TAPE  BEFORE  THE 

C PRESENT  OUTPUT  IS  WRITTEN » WHICH  SHOULD  BE  EQUAL  TO  MAXHR 

C OF  THE  PREVIOUS  RUN.  (OUTPUT  TAPE  UNIT  IS  A) 

C 

c MAXHR  = MAXIMUM  EXPECTED  NUMBER  OF  HOURS  ON  DATA  TAPE  (UNIT  3)  + JCOUNT 

C JCOUNT  = MAXHR  OF  PREVIOUS  RUN 

C FOR  2ND  DATA  CARD—  (ALL  INTEGERS) 

C MONTH  = MONTH  OF  1ST  DATA  SET  IN  COL  1-3 

C NDAY  = DATE  OF  1ST  DATA  SET  IN  COL  4-6 

C NYR  = YEAR  OF  1ST  DATA  SET  IN  COL  7-9 

C MDAY  = NUMBER  OF  DAYS  IN  MONTH  IN  COL  10-12 

C NDATA  = ALWAYS  100  IN  COL  13-15 

C NTIME  = TIME  OF  1ST  DATA  SET  IN  COL  16-19 

C ONE  CYCLE  OF  DATA  CORRESPONDS  TO  1324  CHARACTERS 

C ONE  RECORD  OF  TAPE  CONTAINS  510  CHARACTERS 

C 

c 

c 

U I MENS I ON I DATA ( 15) /060  »0261 > 0262 » 063 1 0264 » 065 » 066 » 0267 f 0270 » 07l > 

1053 > 0255 » 0240 » 0215 r 0253/ 

DIMENSION  I TABLE (l4)/0»l»2»3»4»5»6»7»8»9»5l»52»53»54/»NR 

2(12) * ICHN (100) * VI (100) »YR(12) » NOTES ( 3Q )» REM ( 100 ) 

DIMENSION  REMK (8) » TEMP(IOO) »HF(10) » IFOR(6000)>  FOR(IOO) 

DATA  ( REM ( JH ) * JH=1 » 100 ) /3HREF » 78*1HT t 2H** » 3+3HS0L t 3HWDD » 3HWDV * 

1 8*2HHF» 

2 1 HT  t 2H  t 4HR A I N » 6H  REF  * 6H  UVR  ,2*3H***/ 

DATA ( REMK ( IK ) » IK=1»8)/2HN  »2HNW»2HW  >2HSW»2HS  »2MSL»2HE  »2HnE/ 

INTEGER  SKlPf  FSTCNT 

INTEGER  TPETST 
COMMON  I 1 1 » NT  APE • NCNT 
I 11=510 
NT  APE  = 8 

CALL  CLOCKT(TIMEO) 

JCOUNT  = 0 

READ <5» 7504)  KCARD >SK IP » MAXHR » JCOUNT » TPETST 
WR I TE (6» 7504) KCARD » SKIP* MAXHR » JCOUNT r TPETST 
7504  FORMAT (1017) 

FSTCNT  = JCOUNT  + 1 

-READ (5» 10)  MONTH » NDAY » NYR » MDAY » NDATA »NHOUR» MIN 
WRITE(6» 10) MONTH r NDAY » NYR » MDAY » NDATA »NHOUR» MIN 
10  FORMAT ( 51 3 » 2 12 ) 
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CALL  TPSK IP ( 7 » SK IP ) 

I u'  A Y = NDAY 
7 Q U 1 FORMAT ( 10F7. 0 ) 

IF  (KCARD.EO. 0 ) GO  TO  7000 
PC  77o6  L = l»2 
LL=(L-1)*13+1 
L J-LL+ 1 2 

Rf AD(b» 13) < NOTES (JL) » JL=LL*LJ) 

7 7 oh  rtH ITE (6» 13)  (NOTES! JL) » JL=LL»LJ) 

IL?LOCK  = 0 

7000  DU  7002  K=1#NDATA 
7 0 UD  I CHN ( K ) — K - 1 

IF (KCAPD.NE.O)  GO  TO  9114 
7004  DO  7003  L=l*2 
LL=(L-1 >*13+1 
Lu-L L+12 

REAP (5. 13) ( NOTES ( JL) » JL=LL»LJ) 

7 0 03  to 9 I T E ( 6 » 1 3 ) ( NOTES ( JL )» JL=LL » LJ ) 

13  FORMAT  (1H1»  1 3 A6  » A2 ) 

RE Ap(5»7001) (VI <K)  » K = 1 # NDAT  A ) 

GC  TO  7006 
9114  CONTINUE 

CALL  CLOCK 
r, TESTED 

LSTNO  = LSTNO  + 1 
JCOUNT  = JCOUNT  + 1 
IF ( JCOUNT . GT . MAXHR ) GO  TO  69g 
!*• TEST  = 0 
27  IN-0 

2 < - CmLL  CHAR(N) 

IF ( III. EQ. 510)  IBLOCK=IRLOCK+l 
N TEST  = MTEST  + 1 
752U  FORMAT ( 1012*2110) 

7516  CALL  CLOCKS ( TIME ) 

T I mf=j I ME-T I MEO 
IF ( T IMF -50 • ) 75 1 0 » 75 1 1 » 751 1 
7511  IF (MTEST. LT. 1324)  GO  TO  7510 
7522  toR I TE ( 6 » 7512 ) MTEST 

VvR  I TE  ( 6 » 7520 ) N»III,IBLOCK 

7 j 1 2 FORMAT ( 1 H 1 * T 00  BAD  YOU  ARE  IN  TROUBLE  » TAPE  IS  RAD*  MTtST='I*0> 
GO  TO  31 

7510  IFOR (MTEST )=N 

IF (N.NE. IDATA ( 1 ) ) GO  TO  27 
IN=IN+1 

II-  ( IN.LT.2)  GO  TO  2B 
85t>5  CALL  CHAR(N) 

MTEST  = MTEST  + 1 
IFOR (MTEST)  = N 
I F ( N . EG • IDATA ( 1 1 ) ) GO  TO  7010 
IF (N.EO. IDATA ( 15) ) GO  TO  7010 
IF  (N.EO. IDATA ( 1 ) ) GO  TO  8585 
IF ( IN.LT.3)  GO  TO  28 
GO  TO  27 
7ulL  CONTINUE 
7008  LST*)0=n 
70  u9  DO  4 J-  1 • 9 

CALL  CHAP ( M ) 

DO  3 1=1*12 
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IF (N.NE. IDATA(I ) ) GO  TO  3 
IMR(J)  = ITABLE(I) 

YR ( J)=NR( J) 

GO  TO  4 

3 CONTINUE 

4 CONTINUE 

VI ( 1)=YR(1)*1OOO.+YR(2)*1OO.+YR(3)*1O.+YR(4)+YR(5)/1O+YR(6)/10O. 
7b U i I CHI'. ! ( 1 ) = 0 

42  DO  100  K=2*NDATA 
DO  101  J=l*13 
CALL  CHAR ( N ) 

DO  103  1=1*12 

xr  iUATA  ( I ) ) GO  TO  103 

NR ( J ) =1 T ABLE ( I ) 

YR ( J)=NR( J) 

GO  TO  101 
103  CONTINUE 
101  CONTINUE 

I CHN ( K ) =NR ( 1 ) *1 00+NR ( 2 ) *1 0+NR ( 3 ) 

ICHNX=ICHN(K)-ICHN(K-1) 

VI (K)=YR(5)*10+YR(6)+YR(7)/I0.+YR(8)/I00.+YR(Q)/I000+YR(10)/Iu000 
IF (NR (4) .EQ. 52)  V1(K)=-  VI (K) 

100  CONTINUE 

IF (TPET5T.EQ.1 ) WR I TE ( 6 r 752 1 ) ( Vl C K ) *K  = 1 *NDATA) 

7521  FORMAT ( 1 0F1 0 • 4 ) 

C 

L CHANGE  MILLIVOLTS  TO  DEG  F 

C 

7006  TEMP ( 1 ) =1000  . 

DO  99  I = 2 * 68 
PV  = ARS( VI ( I ) ) 

V2  = 2.711  - PV 
TEMPI I)=TCONF(V2) 

99  CONTINUE 

DO  98  T = 76*79 
98  TEMPI  I ) = TCONFIVl  ID) 

TEMP  I 69)  - TCONFIVl (69) ) 

DO  97  I : 70*75 
V 2 = V 1 ( I ) + VI (69) 

97  TEMPI  I ) = TCONF ( V2 ) 

FOR  (81)  = 99 . 4* VI  I 81 ) 

FOR  I 82 ) = 12  « 6*V1 ( 82 ) 


DETERMINE  HEAT 

FLUX 

HF(l)  = 

1.44 

* 

VI (86) 

HE (2)  = 

1 .40 

* 

VI (87) 

FIFI3)  = 

1.37 

* 

VI 188) 

HF I 4 ) = 

1.44 

* 

VI (89) 

HF I 5 ) - 

1.51 

* 

VI (90) 

HF I 6 ) = 

1.44 

* 

VI (91 ) 

HF  I 7 ) = 

1.37 

♦ 

VI (92) 

HH8)  = 

1.37 

* 

VI (93) 

DETERMINE  WIND 

VELOCITY 

WOV  = 214.  * V 1 ( 85 ) 
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DETERMINE  WIND  DIRECTION 


C 


VvLD  = VI (84) 

IF (WHO. GE. .350)  GO  TO  61 
IK  r l 

00  TO  93 

ol  IF  (WDiJ.GE.  .450.  AND.WDD.LT.  .550) 
I K - ? 

GO  TO  9 3 

oj  IF  ( '•  00  . GE  . .650.AND.WDD.lt.  .770) 
If  =3 

GO  TO  93 

o3  If  (wDU.GE. .900. AND.WDD.LT. 1 . 10 ) 
I K — 4 

GO  TO  °3 

u4  I F ( v;DD  . GE  . I . 50  . AND  . WDD  . L T . 1 . 90  ) 

1 ft  — 5 

GO  TO  93 

03  IF  (i-OJ.GE.  1.15.  AND.  WDO.LT.  I .30) 
I k =6 

Gu  TO  -N3 

Ob  If  ('•  DO.  GE.  4. 50.  AND.  WDD.LT.  6. 00) 
IK  = 7 

GG  TO  Qf 

o7  If  ( '•  DD.GE.2 .50  . AND . WDD.LT  . 4 . 40  ) 


GO  TO  6? 


GO  TO  63 

GO  TO  64 

GO  TO  go 

GO  TO  66 

GO  TO  67 


IKzR 


v SET  TtU-'P#  HF # SOLAR#  WIND  INTO  VARIARLE  CALLED  FO:<  ( INF  OKMi.T  „ C,  ) 

c 

95  LC  400  JB=1»79 
4 u U FOR(Jri)  = TEMP (JR) 

FOR (80)  = VI  (80) 

F OR ( 83 ) = 0. 

FOR ( 84 ) =REMK ( IK ) 

FOR (83>  =WOV 
DO  GOu  JC- 1 » 8 
JF-86+JC 

oUO  FOR(JB)  = HF(JC) 

FOR (94)— TCONF (2.711— ABS(V1 (94) ) ) 

F OR ( 95 ) -VI (95) 

F Of  (96)=V1 (95)*200./Vl (96) /1750. 

PV  r VI (97) 

V 2 = 2.711  - PV 
FOP (97)  = TCONF  (V2) 

FCR(98)=V1 ( 98 ) /0 . 925 
CAl  L CL  OCKS(TINE) 

T I ME- T I ME”T IMEO 
IF ( TIME-150. ) 7513  » 7514  # 7514 
7 3 i 4 WR ITF ( 6 » 751 5 ) TIME 

7515  FORMAT ( 1H0 * SOMETHING  IS  WRONG  WITH  OAT  A # FL  APSFD  TIME  ThIS  KEtC-n 
IIS'  # F 1 0 . 4 ) 

GO  TO  70 
G 

C WRITE  OUT  ANSWERS  IN  TABLE  FORM 

G 

7313  IF  (KCARD.EQ.O)  GO  TO  43 

W R I T F ( 7 ) JCOUNT  »LSTNO»  MONTH  »NDAY  »NYR  fNHOIJR  »M  In#  ( f Ok  ( J ) # o-l  # 1 0 u ) 

WP I TP ( 6 # 7502 ) JCOUNT » LSTNO 

73u2  F C'NA  T ( 8H  RUN  NO=1I10#18H  PRINTER  LIST  NO  = »IilO  ) 
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non  n o r n o n n non 


IF (MHOUR.NE.O)  GO  TO  7503 
GO  TO  70 

7503  IF  (flHOUR.NE.  12)  GO  TO  PI  13 
70  WRITE (6 » 20)  MONTH . IDA Y » NYR 
20  FORMAT  ( ’1*  . ’WEEK  OF  *.313) 

WRITE (6. 21 ) 

cl  FORMAT ( IX » ’REMARKS* ) 

<+5  Wk  I TF  ( 6 » 30 ) MONTH  . NDA  Y » NHOUP  » MlN  ► JCOUNT 

50  FORMAT  (3X*  ’DATE’  . I3»  ’ /’.  13.  4X . * TlMfT  ’ . 2X . 2 1 2 » 4X r » RUN  uO.»» 

1 15) 

WRITE (6» 180)  LSTNO 

lofi  F0RMAT(33X. ’PRINTER  LIST  NO  =’.110/) 

43  WRITE  (6.187) 

187  FORMAT  (IX.  2(’CHAN  NO ' . 5X » • MV  * » 5X » ’ DATA  REMARK’.  4X)) 

DO  1200  IL=1 . 33 

12U0  WRITE (6 .201) ICHN( IL) »V1 ( IL) »FOR( ID » REM ( TL ) . ICHN( IL+bO) » VI  ( Id  + * 0 ) . 
1 FOP( IL+50) »REM( IL+50) 

201  FORMAK  2(I6.F11.4.  F7.1.  ?X.  A5.  4X)) 

WRITE (6.203)  ICHN(34) .VI (34) . FOR(34).  RFM(34).  ICHN(84)»  Vl(64), 

1 FOR (84).  REM (84) 

203  FORMAT ( I6.F11 .4.  F7.1»  2X.A5.  4X.  16.  F11.4.  5X.1A2.  2x.  Ab»4X) 

DO  1201  I L = 35  » 50 

1201  WRITE(6.201)  ICHN ( IL) . VI ( IL ) » FOR(IL).  REM(IL).  1 CHN ( Il+50 ) . 

1 V 1 ( IL+50) . FOR ( IL+50 ) . REM(IL+50) 

WRITE (6.204) 

2u4  FORMAT (1H1) 

IF (TIME -50. ) 7516.7516.7517 
7516  IF (KCARD.EQ.O)  GO  TO  7000 
C ADVANCE  TIME  BY  ONE  HOUR 

C 

9113  NHOlJR  = NHOUR  + l 

IF (NHOUR.LT. 24)  GO  TO  31 


START  NEW  DAY 

nhour=o 

ADVANCE  DATE  BY  ONE 

NLA Y = NDAY  + 1 

CHECK  FOR  LAST  DAY  OF  MONTH 

SET  MO AY  EQUAL  TO  NUMBER  OF  DAYS  IN  MONTH 

IF (NDAY . LE • MDA Y ) GO  TO  31 
NO AY  = 1 

START  NEW  MONTH 

MONTH  = MONTH  + 1 
IF (MONTH. LT. 13)  GO  TO  31 

MONTH  = 1 

nyr=nyr+i 

51  CONTINUE 

MXCNT = JCOUNT -1 
NCNT=MXCNT-FSTCNT+1 
If (TIME-50.)  7523. 7523. 6P9 
7623  IF (KCARD.EQ. 1 ) KPAT=1 
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1MKCARD.EQ.0)  KPAT=? 
GO  TO  (9114,7000)  » K P A T 
699  CALL  CLOCK 
NCNT=NCNT-H 
DO  696  J- 1 * NCNT 
696  BACKSPACE  7 

CALL  TAP(NCNT) 

7617  STOP 
END 
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SUBROUTINE  CHAP(N) 

COMMON  I 1 1 » NT APE  » NONT 
DIMENSION  INCOME (510) » INPUT( 170) 

IF ( II I.NE.510)  GO  TO  2 

NEOF  - 0 

12  CALL  NTR AN ( NT  APE  r 2 » 1 70  » T NPUT » L ) 

1 IF  (L+l ) > 1 
II  (D  5 

1 = 1 

DO  10  J= 1 » 51 0 » 3 

INCOME* J) = ( AND  ( FLD  ( 6 » b » INPUT  ( I ) ) » 1 5 ) *16 ) + 

1 ( AND ( FLO ( 0 » 6 » INPUT ( I ) ) » 1 5 ) ) 

INCOME* J+l ) = ( AND ( FLD ( 18»6> INPUT ( I ) ) r 15) ♦lb)  + 
1 ( AND  ( FLD  ( 12*6*  INPUT  ( I ) ) r 1.5)  ) 

INCOME ( J+2)=(AND( FLO (30 »6» INPUT (T ) ) *15) *16) 

1 + ( AND ( F LD ( 24  t 6 » INPUT ( I ) ) » 1 5 ) ) 

10  1=1+1 
I 11  = 0 

2 IT  1 = 111  + 1 

N= INCOME (III) 

RLTI IRN 


5 ? j = — 1 

WF I T F ( 6 » 1 1 ) L 
I F C L . NF.-2)  GO  TO  13 
IF (NEOF .GT.2)  GO  TO  13 
NL  OF  = NEOF  + 1 
CALL  NTR AN  (NTAPE»22) 

GO  TO  12 

11  FORMAT  (17H0  T jaPE  ERROR  L=I10) 
13  IF (NCNT.EG.O)  STOP 
DO  698  J= 1 » NCNT 
598  BACKSPACE  7 

CALL  TAP (NCNT) 

STOP 

END 
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SUBROUTINE  TAP(NRUN) 

DIMENSION  FOR (100)  » REMK(O) 

DATA  ( RE'MK  ( I K ) * IK=1*8)/2HN  *?HNW*2HW  »?HSW*2H?  *2HSi_»2HE  * 2HNT.  / 
DO  ? >1=  If  NR  UN 

READ ( 7 ) JCOUNTf LSTNOr MONTH rNDAYfNYR f MHO*  «P *wiN*  (TOR ( J) # J-l f 100 ) 

WR l TE ( 6 f 3 ) JCOUNT f LSTNO * MONTH* NDAY , NYP f NHOi  JR  R'p: 

5 FORMAT (1017) 

DO  9 IK=1*Q 

THE  TEST  FOR  EQUALITY  BETWEEN  NON- I NTE^ERS  MAY  NOT  ciL  MEANINGFUL. 
IF (FOR (84) .EO.REMK (IK) ) FOR ( 84 ) -IK 

9 CONTINUE 

2 WRITE (6 f 4 ) (FOR (K ) *K  = 1 * 100 ) 

4 FORMAT (1  OF  10.1) 

re  turn 

END 
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G TT  1'PEKATURF  (IN  DEGREES  CELSIUS)  USING  A FOURTH  DEGREE  EQUATION 

C THIS  SUBROUTINE  WILL  CONVERT  A THERMOCOUPLE  EM r ( in  MILLIVOLTS) 

C FITTED  TO  TABLE  10  IN  MBS  C56 1 ( COPPER-CONSTANT A!  i ) 

C IT  TCOMF  IS  USED  AS  THE  ENTRY  THE  TEMPER A TURF  WILL  bL  IN 

C DEGREES  FAHRENHEIT 

FUNCTION  T CON  (X) 

DIMENSION  A(S)»  B<5) 

DATA  (A(N)  » N — 1 » 5)  / 13. 2 028* -35-5181 » 46. 33851  » 38. 48602-»  0 . / » 

1 (B(M)  rN=l  »5>/0.  >52.81 12  » -1  06 . 5543»  Q2 . 67702  * 38.48-b02/ 

M=0 

i IF  < (X.GT.-5.4) .AND. (X.LT. 20.9) > GO  TO  3 
T C ON  - 9999.9 
Rt TURN 

3 T - X ♦ ( 0. 0314—0 • 000755+X ) 

4 E = A ( 1 ) 

S-R ( 1 ) 

LO  5 J-2fb 
E-A ( J) +T*E 

5 S=ti(J)  + T*S 

n - x - e 

IF  (ABS(D) .LT.l.E-6)  GO  TO  10 
T=T fO/S 

GO  TO  4 

i ) TC ON- 1 0 0 0 . * T 

IF  (M.FQ.l)  TCON=l .8*TCON+32. 

R L TURN 

El  TRY  TCOMF(X) 

M - 1 

GO  TO  1 
r IT 


i N T 0 
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SUBROUTINE  TPSK IP ( UNIT  » SK I P ) 
INTEGER  UNIT  # sk  ip 
II  (SK IP) lr 2r 3 
i SKIP  - -SKIP 
DO  4 J- 1 » SKIP 

4 BACKSPACE  UNIT 
GL  TO  2 

o DC  5 J=1 > SKIP 

5 READ (UNIT) 

2 Rt TURN 

FDD 
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Attentions  Mr.  D.  M.  Trwnble 


Subject:  Contract  H67274  (S-442247-68) 

Subtoil  Investigation  and  Soil  Tatting, 
Earth  Tesperature  Iattruawntation 
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GENERAL  NOTES 


1.  Number  in  right  hand  column  of  boring  log  indicates  blows  required  to  drive  ?r.  In.  O.D.,  . In.  1.0. 

sampling  spoon  one  foot,  using pound  hammer  falling  3 inches. 

2.  Classification  of  soil  is  by  visual  inspection  and  is  in  accordance  with  the  Unified  Soil  Classification  System. 

3.  Estimated  ground-wafer  levels  indicated  by  Y ; these  levels  are  only  estimates  from  available 

data  and  may  vary  with  precipitation,  porosity  of  the  soil,  site  topography,  etc. 

4.  BOPAMG  OWUJtO  fe\  MECUNMICAL.  NOSER  \MlTUOUT  USE  Of  ORIUAUG  VvJKTER. 

5.  B0R\MG  LKVOUT  FUR.MISHEO. 
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